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Abstract. Magnetic fields play a crucial role at all stages of the formation of low 
mass stars and planetary systems. In the final stages, in particular, they control the 
kinematics of in-falling gas from circumstellar discs, and the launching and collimation 
of spectacular outflows. The magnetic coupling with the disc is thought to influence the 
rotational evolution of the star, while magnetised stellar winds control the braking of 
more evolved stars and may influence the migration of planets. Magnetic reconnection 
events trigger energetic flares which irradiate circumstellar discs with high energy 
particles that influence the disc chemistry and set the initial conditions for planet 
formation. However, it is only in the past few years that the current generation 
of optical spectropolarimeters have allowed the magnetic fields of forming solar-like 
stars to be probed in unprecedented detail. In order to do justice to the recent 
extensive observational programs new theoretical models are being developed that 
incorporate magnetic fields with an observed degree of complexity. In this review 
we draw together disparate results from the classical electromagnctism, molecular 
physics/chemistry, and the geophysics literature, and demonstrate how they can be 
adapted to construct models of the large scale magnetospheres of stars and planets. 
We conclude by examining how the incorporation of multipolar magnetic fields into 
new theoretical models will drive future progress in the field through the elucidation 
of several observational conundrums. 
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1. Introduction 

The current generation of spectropolarimeters, ESPaDOnS at the Canada-France- 
Hawai'i telescope, and its twin instrument NARVAL at Telescope Bernard Lyot in 
the French Pyrenees, are revolutionising our understanding of stellar magnetism as a 
function of stellar age and spectral type. Results include (but are not limited to) the 
possible detection of a remnant fossil field on a hot massive star |49j; the first ever 
magnetic surface maps of pre-main sequence stars in the classical T Tauri phase of 
their evolution [50] EH EO H07j ; the discovery of successive global magnetic polarity 
switches on a star other than the Sun, whose short cycle may be caused by the known 
presence of an orbiting close-in giant planet [53] ; the rapid increase in field complexity 
at the transition from completely convective low-mass stars to those with radiative cores 
[1801 151] ; and the discovery of globally structured magnetic fields on the intermediate 
mass Herbig Ae-Be stars [29] . 

Knowledge of the medium and large-scale topology of stellar magnetospheres 
provided by the spectropolarimetric observations is crucial to our understanding of 
many important processes. For low-mass pre-main sequence stars the magnetic star- 
disc interaction is believed to control the spin evolution of the central star, and may 
also be responsible for the collimation and launching of outflows [1701 H73|. [66] 1147] . 
As both low and high mass stars evolve on the main sequence, the angular momentum 
that can be extracted by winds depends on the amount and distribution of open field 
[103;, 1249] . Orbiting close-in giant planets may also interact magnetically with stellar 
magnetospheres, which in principle provides a mechanism for characterising planetary 
magnetic fields and therefore their internal structure [183^ I227| I114j . 

In order to model such physical processes, new theories and simulations that 
incorporate magnetic fields with an observed degree of complexity are required. Over 
the past few years a series of models which move beyond the assumption that stellar 
magnetic fields are dipolar have been developed. In this paper we review such models, 
and provide a thorough derivation of the magnetic potential in the region exterior to a 
star, deriving general expressions for a large-scale multipolar stellar magnetosphere. In 
this review we concentrate on the theoretical study of stellar magnetospheres, briefly 
discussing observational results where appropriate. An authoritative review of the 
observational study of stellar magnetic fields is provided in [56]. We focus on the 
magnetic fields of forming solar-like stars, although the analytic expressions derived 
herein are applicable to models of stellar and planetary magnetospheres generally. In 
the remainder of §[T] we review the basic techniques that allow stellar magnetic fields to 
be detected and mapped, before focussing specifically on the magnetic fields of accreting 
T Tauri stars - low mass stars still surrounded by planet forming discs. In ^2] we discuss 
the observational support for the magnetospheric accretion scenario and briefly review 
previous models with dipolar stellar magnetic fields. Following this, in £j3j we draw 
together results from molecular physics and classical electromagnetism to derive self- 
consistent analytic expressions for multipolar stellar magnetic fields. In ^4] we discuss 
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the first models of the magnetospheres of forming solar-like stars that take account of 
non-dipolar magnetic fields. We conclude in §[5] by highlighting several open problems 
where consideration of the true complexity of stellar magnetic fields may be crucial for 
future progress. 

1.1. Detecting and mapping stellar magnetic fields 

Stellar magnetic fields can be probed using two complementary techniques. Measuring 
the Zeeman broadening of unpolarised spectral lines has proved to be a successful 
method of determining average stellar field strengths. References 120 lj and |205] 
demonstrated that by measuring changes in the shapes of magnetically sensitive lines in 
intensity spectra, estimates of total field strength, and the fraction of the stellar surface 
covered in fields (the magnetic filling factor) could be estimated. Zeeman broadening 
measurements were carried out on a number of stars (for example [2144 1118] ). however, 
for young stars this proved problematic due to rotational broadening dominating the 
line profiles [12]. Broadening measurements are easier to carry out at infrared (IR) 
wavelengths, as the Zeeman splitting increases more rapidly at longer wavelengths 
compared to Doppler broadening [121] . The use of IR line profiles to measure stellar 
magnetic fields was pioneered by the authors of |216] and |215] . Subsequently the 
analysis of various features in IR spectra has proved to be an extremely successful 
method of detecting stellar magnetic fields (for example [20 0], I201[ [124J). Zeeman 
broadening measurements, however, give no information on the magnetic field topology. 
In contrast, measuring the circular polarisation signature in spectral lines gives access to 
the field topology (see [3B] and [SO] for reviews of the basic methodology). However, like 
all polarisation techniques, this suffers from flux cancellation effects and yields limited 
information regarding the field strength. 

If a stellar atmosphere is permeated by a magnetic field, spectral lines forming 
in that region will be polarised, with the sense of polarisation depending on the field 
polarity. In practise the polarisation signals detected in photospheric absorption lines 
are small, and cross-correlation techniques (such as Least-Squares Deconvolution; [55] ) 
are employed in order to extract information from as many spectral lines as possible. 
The signal-to-noise ratio of the resulting average Zeeman signature is several tens of 
times larger than that of a single spectral line [55]. Magnetic surface features produce 
distortions in the Zeeman signature that depend on the latitude and longitude of 
the magnetic region, and on the orientation of the field lines. By monitoring how 
such distortions move through the Zeeman signature as the star rotates, a method 
referred to as Zeeman-Doppler imaging (ZDI; |226j), the 2D distribution of magnetic 
polarities across the surface of stars can be determined using maximum entropy image 
reconstruction techniques [23J. The field orientation can also be inferred within the 
magnetic regions [33]. After the initial success of ZDI in recovering the first magnetic 
maps of a star other than the Sun [30], the technique has been applied to construct 
magnetic surface maps of stars of different spectral type at various evolutionary stages 
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(for example, pEEl Ell SSI ED E3 ES EED]). In the latest incarnation of ZDI the field 
topology at the stellar surface is expressed as a spherical harmonic decomposition |49j. 
The surface field is described as the sum of a poloidal plus a toroidal component, which 
allows us to determine how much the field departs from a purely potential (poloidal) 
state. For ZDI of accreting T Tauri stars polarisation signals in both accretion related 
emission lines and photospheric absorption lines are considered when constructing 
magnetic maps [SUJ [Ml EI]- Photospheric absorption lines form across the entire star, 
while accretion related emission lines form locally where accreting gas impacts the stellar 
surface. Thus magnetic maps constructed from photospheric absorption lines only miss 
strong field regions which contain the foot points of the large scale field lines that 
interact with the disc [50]. The resolution of such magnetic surface maps is limited by 
the stellar rotation period and inclination, with a finer spatial resolution at the surface 
of the star achievable for faster rotators, and by the flux cancellation effect of two nearby 
opposite polarity regions giving rise to oppositely polarised signals, resulting in a zero 
net polarisation signal (252]. As a result, the smallest scale magnetic features, such as 
stellar equivalents of the small bipolar groups observed on the Sun (for example, [4T]). 
remain unresolved. Instead, spectropolarimetric Stokes V (circular polarisation) studies 
are limited to probing the medium and larger scale fields, and likely miss a large fraction 
of the total magnetic flux [2U2j . None-the-less, significant advances in the study of stellar 
magnetism have been made over the past few years using optical spectropolarimeters, 
in particular in the mapping of the magnetic fields of forming solar-like stars, as we 
discuss in the following subsection. 

Zeeman-Doppler imaging is not the only method that has been developed to map 
stellar surface fields. Magnetic Doppler imaging (MDI), which can trace its origins back 
to work that pre-dates the development of ZDI [194L 178] . is an alternative technique 
that incorporates polarised radiative transfer, and can also include linear polarisation 
diagnostics (Stokes Q and U) in the field reconstruction |195[ I131[ I133[ I196[ 1132] . 
However, MDI has thus far only been applied to construct maps of a few chemically 
peculiar stars |13U[ llb'5[ I134j . Furthermore, as argued by the authors of [51], MDI 
is currently limited to only the brightest and most magnetic stars. None-the-less, 
development of MDI will continue to be scientifically fruitful in future years, and will 
provide important comparison tests with the results of ZDI studies. 

1.2. Accreting T Tauri stars and observations of their magnetic fields 

Classical T Tauri stars (CTTS) represent a key transitional period in the life of a star, 
between the embedded protostellar phase of spherical accretion and the main sequence 
stage. They are low mass pre-main sequence stars which accrete material from dusty 
circumstellar discs. They possess strong magnetic fields, of order a few kG |124] , which 
truncate the disc and force in-falling gas to flow along the field lines. Material rains 
down on to the stellar surface, where it shocks and produces hot spots that emit in the 
optical, ultraviolet (UV), and X-ray wavebands. CTTS are observed typically to rotate 
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well below break-up speed, and are more slowly rotating (on average) than older pre- 
main sequence weak-line T Tauri stars whose discs have largely dispersed (for example, 
[T7]). CTTS can be in excess of 1000 times more active in X-rays than the Sun is 
presently. X-rays from the central star may influence the dynamics and chemistry of 
the circumstellar disc, which will in turn set the initial conditions for planet formation 
[63, 190J. Understanding the final stage of formation of CTTS, and how they interact 
with their discs, is crucial if we wish to understand the formation of the Sun and our 
own Solar System. Many accreting T Tauri stars will eventually evolve along the main 
sequence surrounded by planetary systems much like our own. 

There is an abundance of observational evidence which supports the basic scenario 
of magnetically controlled accretion from truncated circumstellar discs. Excess IR 
emission is consistent with CTTS being surrounded by dusty discs, while the shapes 
of spectral energy distributions (SEDs) in the near-infrared (nIR) are consistent with 



detected in many emission lines, with broad redshifted absorption components indicating 
gas infall at approximately free-fall velocities, consistent with columns of gas being 
magnetically channelled on to the stellar surface [HUES]- Blue shifted absorption is also 
commonly detected, indicating that strong outflows are common from CTTS, although 
whether such outflows originate mainly from the star, or from the disc surface remains 
an open question [147] . Excess continuum emission at optical and in particular UV 
wavelengths is consistent with the existence of accretion shocks at the stellar surface, 
formed due to the high velocity impact of accreting gas. This excess emission, the veiling 
continuum, makes absorption lines shallower that they would appear in non-accreting 
stars of the same spectral type |99J . Estimates of the amount of veiling provides a method 
of determining the mass accretion rate on to the star [9?J. CTTS are highly variable at 
all wavelengths, over a variety of timescales. Such variability is thought to arise due to 
a complex mixture of hot (accretion related) and cool (magnetic flux emergence related) 
spots distributed across the stellar surface, time variable mass accretion and outflows, 
as well as obscuration effects such as warped inner discs, and columns of accreting gas 
rotating across the line-of-sight [18]. Meanwhile copious amounts of X-ray emission, 
thought to arise due to particle acceleration along field lines following reconnection 
events, indicates that CTTS are extremely magnetically active (see [H] and |177] for 
comprehensive reviews) . 

The magnetospheric accretion scenario requires that T Tauri stars possess strong 
magnetic fields that are sufficiently globally ordered to truncate the disc. Measuring 
their magnetic fields, however, remains difficult as the stars are faint, and subject to high 

| Inner dust disc radii derived from interferometric measurements have often been found to be larger 
than that derived from SED fitting, see j!66j and references therein. This discrepancy may arise from 
the crudeness of the disc models used to convert from interferometric visibilities to disc inner radii 
[193] . For CTTS, however, it is the location of the inner gas disc, which extends beyond the dust disc, 
for example pQ, that is important. Interferometric studies are just beginning to probe gas on such a 
small scale, for example [61] (see also the CO transition spectroscopy work in |185j and [28] ). 




Inverse P-Cygni profiles are commonly 
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levels of spectral variability. Initial spectropolarimetric studies at optical wavelengths 
failed at directly detecting magnetic fields (221 H26[ 1127] . however, due to the flux 
cancellation effect whereby signal from regions of opposite polarity cancel, these initial 
failures were in fact early indications of the complex nature of T Tauri magnetic fields. 
The first field detections were made by estimating the increase in line equivalent width 
that arises due to the saturation of the Zeeman components. Reference [96], through 
careful analysis of photospheric Fe I absorption lines, found the product of magnetic field 
strength and (magnetic) filling factor of order ~kG on a few accreting T Tauri stars 
(reference [12] had previously used the same technique to make similar field detections 
on non-accreting T Tauri stars). The most successful method of measuring T Tauri field 
strengths has been through the analysis of magnetically sensitive lines at IR wavelengths, 
as Zeeman broadening increases more rapidly with increasing wavelength compared to 
rotational broadening. Strong average fields of order 1-3 kG have now been detected 
on T Tauri stars in Taurus, the Orion Nebula Cluster (ONC) and the TW Hydrae 
Association p2H EMI EEH1 [263 E52J [121]. Typically magnetically sensitive IR lines of 
Ti I are used, their shape being best described if the stellar surface contains a distribution 
of field strengths (up to ~ 6 — 7kG). Averaging over the distribution yields photospheric 
surface fields of ~ 1 — 3kG, with similar values for both accreting and non-accreting (i.e. 
disc-less) T Tauri stars |123l 1252] . Average field strengths of this magnitude are strong 
enough to disrupt the disc (as discussed in the following section), however, reference 
[217] points out that such strong fields are not necessarily sufficiently globally ordered 
to truncate the disc at several stellar radii. 

Zeeman broadening measurements of lines in intensity spectra have yielded several 
intriguing results. The average field strengths of T Tauri stars are a few kG. The authors 
of |124j and |267] argue that such strong fields may inhibit the coronal X-ray emission. 
The quiescent X-ray emission is thought to be due to many small flares triggered by 
by reconnection events arising from the release of magnetic stresses built up due to the 
motion of field line foot points. Strong fields may inhibit the foot point motions and 
the consequent tangling of the field. This may explain why T Tauri stars appear to 
be less luminous in X-rays than predicted from a correlation between X-ray luminosity 
and (unsigned) magnetic flux 11241 1266] . Another major finding from Zeeman 
broadening studies of stars in different star forming regions is the apparent decrease in 
unsigned magnetic flux (4iiRlB) with age [267] . Such a trend remains unexplained. 

T Tauri stars, the majority of which are completely convective, host magnetic fields 
that are most likely to be dynamo generated, see for example |221[ 1143] , and [131 [21] , 
and references therein, for some of the recent work on magnetic field generation in 
completely convective stars in general. However, it is occasionally suggested that their 
magnetic fields may possess a fossil component |245l 159"] I182[ IT 24] . Fossil magnetic fields 
are fields that have survived from the initial collapse of the magnetised cloud core during 
the earliest stages of the formation of the star. The arguments against the existence of 
fossil fields at typical T Tauri ages (~ few Myr) have been succinctly summarised by the 
authors of [56] . Firstly, the onset of convection is thought to rapidly destroy any remnant 
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fossil field on a timescale of not more than 1000 years [31]. Secondly, indicators such 
as flares (commonly observed on T Tauri stars) suggests reordering of their magnetic 
fields, and thus they cannot be linked to evolutionary processes from millions of years 
in the past. Thirdly, the similarity between the large scale magnetospheric properties 
of T Tauri stars and those of M-dwarfs (with an age of order Gyr i.e. so old that 
their fields are certainly not fossil), that we discuss at the end of this section, is further 
evidence that fields are dynamo generated. Although dynamo magnetic field generation 
models for partially and fully convective stars are still debated, the current generation of 
spectropolarimeters is providing the community with crucial information on how stellar 
field topologies vary with quantities such as stellar mass and rotation period (see the 
review article [56]). 

Zeeman broadening measurements do not yield information about the magnetic 
geometry of accreting T Tauri stars. However, small wavelength shifts in spectral lines 
observed in right and left circularly polarised light provide another method of diagnosing 
stellar magnetic fields. As previously mentioned, the earliest spectropolarimetric studies 
failed at detecting T Tauri magnetic fields. A major break through was the detection 
of strong circular polarisation in the He I 5876A line [120J. This line of He I has a 
high excitation potential and is thus believed to form at the accretion shock, where 
columns of magnetically channelled gas impact the stellar surface [T4]. Polarisation 
detections in this line are thus tracing the field on the star where the large scale 
field lines that interact with the disc are anchored. The polarisation signal is often 
found to be rotationally modulated, with the variation in the derived line-of-sight (or 
longitudinal) field component with rotation phase well fitted by a simple model with 
a single accreting spot in the visible hemisphere |252j . This is fully consistent with 
the findings from ZDI studies, discussed below, where evidence for single dominant 
accretion spots at high latitudes is found. However, despite arguments that variations 
in the longitudinal field component, derived from polarisation detections in the He 
I 5876A line, were attributable to rotational modulation (for example, |251|, 1252] ). 
[32] refutes such suggestions and argues that the field in the line formation region 
is constantly evolving and restructuring on a timescale of only a few hours. The 
ESPaDOnS/NARVAL spectropolarimetric data presented in [54J, however, clearly show 
that although the He I 5876A line is subject to intrinsic variability, its temporal evolution 
is dominated by rotational modulation. This suggests the magnetic field in the He 
I line formation region remains stable on timescales of longer than a rotation cycle, 
and that T Tauri magnetospheres remain stable, at least over a few rotation periods, 
consistent with earlier line profile variability studies of individual stars (for example, 
[77] I117[ 11191 US]). Strong polarisation signals in He I 5876 A and other accretion 
related emission lines have now been reported on a number of accreting T Tauri stars 
[2511 12521 12351 12431 12651 1501 1541 1571 152] . However, polarisation measurements made using 
magnetically sensitive photospheric absorption lines, which presumably form uniformly 
across the entire stellar surface, yield small longitudinal field strengths, well below 
the average fields obtained from Zeeman broadening measurements |234[ 1235] 1236} 142] . 
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Commonly the net polarisation signal is zero [252J. This suggests that accreting T Tauri 
stars host complex surface magnetic fields. In contrast, the strong (and rotationally 
modulated) polarisation signal detected in accretion related emission lines suggests that 
the bulk (although perhaps not all) of gas accreting on to stars from their discs, lands 
on single polarity regions of the stellar surface. However, it is only in the past three 
years that the geometry of T Tauri magnetic fields have been revealed. 

ZDI studies, combined with tomographic imaging techniques, have now revealed 
the true complexity of the magnetic fields of accreting T Tauri stars. At the time of 
writing surface magnetic maps of six stars have been published, namely V2129 Oph 
[50], BP Tau [M], CR Cha, CV Cha [M], V2247 Oph [57] and AA Tau [52]. All 
have been found to have magnetic fields consisting of many high order components. 
At 1.35 M V2129 Oph is believed to have already developed a small radiative core, 
despite its young age (where the stellar mass has been derived using the Siess et al 
pre- main sequence evolutionary models [233], as with the other stars discussed below). 
The magnetic energy was found to concentrate mainly in a strong octupole component 
of polar strength 1.2 kG tilted by ~ 20° relative to the stellar rotation axis. The dipole 
component was found to be weak, only 0.35 kG at the pole and tilted by ~ 30° relative 
to the stellar rotation axis, but in a different plane from the octupole component. The 
surface field in the visible hemisphere was dominated by a 2 kG positive radial field spot 
at high latitude, with the footpoints of the accretion funnel rooted in the same region, 
but differs significantly from a dipole [50J . Evidence for high latitude (polar) cool spots 
|119] and for high latitude accretion hot spots [242] had already been found previously 
through Doppler imaging of other CTTS. The lower mass and completely convective 
star BP Tau (0.7 M ) is found to have a much simpler field topology with the magnetic 
energy shared between strong dipole (of polar strength 1.2 kG) and octupole (1.6 kG at 
the pole) field components [54] . Two surface magnetic maps were derived for BP Tau 
from circularly polarised spectra taken about 10 months apart, but little variation in 
the large scale field topology was detected. A similar result was found for AA Tau, from 
magnetic maps derived from data taken about one year apart [22] • AA Tau is of similar 
mass, radius, and rotation rate as BP Tau, although its magnetic field is even simpler, 
consisting of strong (~2-3kG) dipole, almost anti-parallel with respect to the angular 
momentum vector of the star, with an octupole component five times weaker. 

The initial ZDI results suggest that the field complexity is intimately related 
to the depth of the convection zone, with completely convective pre-main sequence 
stars hosting simpler dominantly poloidal large scale magnetic fields with strong dipole 
components. These results are consistent with the ZDI study of the massive accreting 
T Tauri stars CR Cha (1.9 M ) and CV Cha (2M ) presented by [107J, as both stars 
have significant radiative cores and have particularly complex large scale field topologies. 
It is also mirrors the findings from ZDI studies of low-mass main sequence M-dwarfs 
[180, 51J. M-dwarfs which are completely convective (those below ~ 0.35 M [30]) were 
found to host simple dominantly poloidal large scale magnetic fields, with strong dipole 
components [180] (the exception being stars below ~ 0.2 M ; only some of which host 
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such simple large scale fields, see below and |181j ). In contrast to the findings for mid 
M-dwarfs, early M-dwarfs (spectral types M0-M3) with small radiative cores have more 
complex large scale fields with strong toroidal and weak dipole components [5T] . Of 
course, due to the flux cancellation that effects circular polarisation studies, discussed 
in the previous section, the dominant and strongest magnetic field components are likely 
be the highest order multipole components that constitute the very small scale field close 
to the stellar surface. This suggestion is emphasised by the authors of [202J who argue 
that the bulk of the total magnetic flux is missed by polarisation studies, with Zeeman 
broadening measurements indicating the presence of small scale field components far 
stronger than those detected by Stokes V studies alone. The work of [192J is also 
consistent with the bulk of the magnetic energy being stored in the small scale features 
that remain unresolved in stellar magnetic maps. 

Zeeman signals are also suppressed within dark (cool) surface spots due to the 
low surface brightness. Cool spots, which on T Tauri stars are believed cover a far 
larger fraction of the stellar surface when compared to sun spots, for example [50] , thus 
represent a potential source of additional missing flux in stellar magnetic maps |125j . 
The flux at the stellar surface is the result of several different processes: the dynamo 
that generated the flux to begin with, the processes that took place during the buoyant 
rise of that flux through the convective zone (and its interaction with the convective 
cells) and finally the surface effects as the flux emerges into the low-density region of 
the photosphere. T Tauri stars are known to have average field strengths of a few kG 
[124] , in excess of the mean solar field strength, although in sun spots the field can reach 
~ 3 — 4kG. It may be the case that the fields in T Tauri cool spots are similarly large 
compared to the mean photospheric field strengths. It is therefore interesting to ask the 
questions: i) what is the relative contribution to the total magnetic flux through the 
surface of T Tauri stars from the dark spots, the flux that is resolved in the ZDI maps, 
and the unresolved flux missing from the ZDI maps due to the flux cancellation effect? 
and ii) are these contributions in the same ratios as we see on the Sun? Untangling the 
different contributions to the total flux through the surface of T Tauri stars promises 
to be a challenge for future theories. 

The picture of completely convective T Tauri stars having simple dominantly 
poloidal large scale fields with strong dipole components may not be valid for the lowest 
mass T Tauri stars. For low mass accreting T Tauri stars (below 0.5 M ) the picture 
may be more complicated. Recently the authors of [57] have presented magnetic maps 
of the low mass T Tauri star V2247 Oph (0.35 M Q ), which has a faster rotation rate 
(Prot — 3.5 d) in comparison with the more massive stars BP Tau (P ro t — 7.6 d) and 
V2129 Oph (Prot = 6.53 d). Various accretion related emission lines are detected in the 
optical spectra, indicating that mass accretion is ongoing in this system, despite little 
IR excess evident from Spitzer satellite data (indicating that the dust component of the 
disc, but not necessarily the gas component, has larger dispersed). The magnetic field 
of V2247 Oph is found to be particularly complex with a very weak dipole component 
compared to that of BP Tau. However, this also appears to be consistent with new 
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results for late-type M-dwarfs (below 0.2 M , or spectral types M5-M8), where individual 
stars are found to host a mixture of complex non-axisymmetric magnetic fields that are 
very different from the strong and simple large scale fields of mid M-dwarfs, and strong 
axisymmetric dipoles which are more consistent with the large scale topologies of mid 
M-dwarfs |181] . The V2247 Oph results demonstrate that more spectropolarimetric data 
for a larger sample of stars are required in order to disentangle the effects of differing 
stellar masses, rotation periods, and mass accretion rates. What is clear, however, is that 
the magnetic fields of accreting T Tauri stars can be significantly more complex than 
a dipole. Before considering how analytic models of multipolar stellar magnetospheres 
can be constructed, we briefly overview the development of T Tauri magnetospheric 
accretion models with dipole magnetic fields. 

2. Magnetospheric accretion models with dipolar magnetic fields 

Although it had been suggested by various authors that T Tauri magnetospheres would 
disrupt circumstellar discs and channel columns of gas on to the star (for example, 
|247l 12481 E7j), it was the inspirational paper of Konigl |140j that demonstrated that 
a multitude of observational features could be explained through the magnetospheric 
accretion scenario. By adapting the Ghosh and Lamb model of accretion on to neutron 
stars [HI [751 ES], Konigl argued that provided T Tauri stars had magnetic fields of 
order a kG that could effectively couple to the disc, discs could be disrupted at several 
stellar radii, the alignment of accretion columns with the line-of-sight could explain the 
development of inverse P-Cygni profiles, while shocks at the base of funnels of accreting 
gas could naturally explain the observed UV excess. Konigl argued that the observed 
slow rotation of accreting T Tauri stars could be explained provided that the spin-up 
torque exerted on the star due to the accretion of high angular momentum material, 
and the magnetic connection to regions of the disc rotating more quickly than the star, 
was exactly balanced by a spin-down torque transmitted by the field lines threading the 
disc exterior to corotation. The corotation radius R co is an important point for models 
of magnetospheric accretion. In the stellar equatorial plane, 



which is the radius at which the Keplerian rotation rate of the disc material is equal to 
that of the star (o>* = 2n/P rot ). Interior (exterior) to this radius, the disc material is 
spinning faster (slower) than the star. At radii interior to corotation, material would 
naturally accrete on to the star, while stellar field lines threading the disc at corotation 
would rotate as a solid body and would not be stretched due to differential rotation. 

Magnetospheric accretion models, such as those proposed in |140| 137] 1188] , provide 
magnetic links between the star and regions of the disc beyond R co which are rotating 
more slowly than the star. By having field lines threading the disc at a range of radii 
the star is able to accrete material without experiencing a net spin-up torque, which 
would act to slowly increase the stellar rotation rate. However, field lines threading 




(2.1) 
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the disc beyond the corotation radius may quickly become wrapped up, inflate, and 
be torn open after only a few rotation periods (see for example jU 12531 D-63J). The 
Shu X-wind model, developed through a series of papers, gets around this problem 
[2291 12301 123T1 12321 ITM H57] . The model introduces the idea of trapped flux, where 
the closed field lines connecting the star and the disc are pinched together in a small 
interaction region about the corotation radius (called the X-region). In such a way the 
strong dipolar field of the star rotates as a solid body with material from the inner part of 
the X-region accreting onto the star. Torques in the funnel flow deposit excess angular 
momentum into the X-region which is then removed by a wind that carries material 
away from the outer portion of the X-region. References |154j and |155j consider the 
funnel flow of gas on to the star, and assuming that accretion occurs at a steady rate, 
also find that the matter angular momentum in the funnel flow is transferred to the 
disc, not to the star, in agreement with the X-wind model. Thus, the Shu X-wind 
model allows accretion to occur without spinning up the star, whilst also providing a 
connection between the accretion process and outflows. Models which combine accretion 
and outflows often predict a correlation between the mass accretion and mass outflow 
rates (for example |187t lB~S]). which has been observed (for example [8]). However, there 
is no reason to expect that discs will always be truncated at the corotation radius. Based 
on IR spectroscopy of CO transitions references |185j and [2H] conclude that gas in the 
inner disc extends to well within the corotation radius, suggesting that there is nothing 
special about corotation (in terms of the location of the disc truncation radius) |§] The 
process of accretion, considered alone, should thus act to spin-up the star, in the absence 
of an additional angular momentum loss mechanism, or significant magnetic connections 
to regions of the disc beyond corotation. 

The author of [259J argues that it is physically impossible in a steady state scenario 
for angular momentum to be transported from the funnel flow region back to the 
outer disc, and thus the material torque must be transferred to the star. Disc-locking 
models (where the stellar rotation rate matches the Keplerian rotation rate at the 
disc truncation radius) have been criticised as being physically [168J unfounded, and is 
often observationally controversial [2381 1239] (we note, however, that other observational 
studies do find good evidence linking the presence of discs, and/or accretion, with slow 
stellar rotation, for example [33J, [62] ) . The authors of [169] argue that for typical T Tauri 
accretion rates, the large scale magnetic field threading the disc is opened to such an 
extent that the star will receive no spin-down torque at all (see also |272j who consider 
the effects of time varying stellar magnetic fields). A recent T Tauri spin-evolution 
model, the first to combine the opening of the large scale magnetosphere due to the 
interaction with the disc with variations in stellar radius and mass accretion rate with 
time, find that all stars experience a net spin-up torque |167j . For their preferred case 
of strong disc coupling, stars are spun-up and end up with rotation periods of less than 
3 days by the end of their simulations at 3 Myr. However, their model only accounts for 

§ The author of [55], who reviews recent progress with the X-wind model, refutes this by pointing out 
that part of the CO emission may come from the accretion funnel itself. 
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spin-down torques provided by a small (in terms of radial extent) connection to the disc 
outside of corotation (with spin-up torques arising from the small connection to disc 
interior to corotation, and from matter accreting on to the star); it does not consider 
additional spin-down torques arising from disc winds, or from stellar winds, that appear 
to be required in order to explain the observed spread in T Tauri rotation periods |167] . 
Clearly there remain many unanswered questions regarding the balance of torques in 
the star-disc system. 

From an observational perspective, reference |100j argues that stellar spin-down 
cannot occur faster than the rate at which angular momentum can be removed by a 
disc wind, or through viscous processes. For stars in the youngest star forming regions, 
such as the ONC studied by the authors of [238], the angular momentum loss time 
scale is comparable to the age of the stars, and thus disc braking may not have had 
sufficient time to slow the stellar rotation rates. Whether-or-not stars can be efficiently 
braked depends crucially on how well the stellar field couples to the disc (for example 
[TOl 11691 fl3l 1167] ). Furthermore, as argued by the authors of [65] (see their section 
2.1), stars must undergo strong braking during the initial optically embedded phase of 
evolution. They [22] argue that the Shu X-wind cannot explain such strong and efficient 
braking, and is thus unable to explain the observed slow rotation of accreting T Tauri 
stars. An alternative magnetospheric accretion/outflow model is the reconnection X- 
wind of Ferreira et al [651 166] where angular momentum is removed by a wind launched 
from the entire surface of the disc. A unique feature of this model is that the outflow 
is powered by the rotational energy of the star itself, and thus the reconnection X-wind 
provides a torque that can efficiently brake the star (see the review in [6?J). Another 
alternative is the accretion powered stellar wind model developed by Matt and Pudritz 
[1701 11721 1173] , which assumes that the spin-up torque due to accretion is balanced by 
the spin-down torque from the stellar wind. However, it is not yet clear how accretion 
can power a stellar outflow (see the discussion in |171j ). A possible suggestion is that gas 
accreting on to the stellar surface in accretion columns excites magnetohydrodynamic 
(MHD) surface waves which drive the stellar outflow [2H]- Unfortunately, the derived 
mass lose rates are an order of magnitude below what is required to extract enough 
angular momentum in the wind to explain the observed slow rotation of accreting 
T Tauri stars. Ultimately the angular momentum removal mechanism may be some 
combination of disc, and accretion powered stellar, winds, for example [39], although 
this is an open question. 

The star-disc interaction, and the process of accretion and outflows, is likely to be 
highly time dependent. In order to incorporate time dependent effects, MHD simulations 
are required. To date, a myriad of MHD models have been constructed of both funnel 
flows and the star-disc interaction, which vary in their assumptions regarding the disc 

physics pm EH EE1 [T351 Effil EEO [HJ [111 IMS [2551 EMU [T59J EST] EEHl E51 1269] . 

Some models predict episodic accretion, periodic inflation, opening, and reconnection 
of the magnetosphere, plasmoid ejection, the launching of winds from the disc, field line 
collimation into jet-like structures, and variable epochs of stellar spin-up and spin-down. 
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The next generation of spectropolarimeters, which will be able to detect magnetic fields 
in the inner disc, will provide data to discriminate between the various MHD simulations. 
A common feature of such models is the assumption that the star possesses a simple 
dipolar magnetic field. The recent ZDI studies discussed in the previous section have now 
demonstrated conclusively that T Tauri magnetic fields are multipolar, with complex 
surface field regions distorting the structure of the large scale accreting field in the 
regions close to the star |94j . 

The authors of [122] took the first steps towards investigating accretion models with 
non-dipolar magnetic fields. They demonstrated that if the dipole field assumption 
is removed from the Shu X-wind model, and under the assumption that the field 
strength does not vary from star to star, there should be a correlation between the 
stellar and accretion parameters of the form Rlf aC c (M*MP rot ) l l 2 (where M is the 
mass accretion rate and f acc the accretion filling factor). Such a correlation agrees 
reasonably well with the observational data. Over the past few years the first models 
of magnetospheric accretion that consider multipolar T Tauri magnetic fields have been 
developed [901 ED E3j M, UM USD HH2l H79l [211]. Before reviewing such models in Q 
we explore results developed for multipole field expansions in classical electromagnetism 
and molecular physics in order to demonstrate how simple analytic models of multipolar 
stellar magnetic fields can be constructed. 

3. Multipole magnetic fields 

Multipole expansions for the potential of a finite static charge distribution in 
electrostatics, and of a continuous current distribution in magnetostatics, are commonly 
encountered in the electromagnetism literature (for example [2S1 E7J 1199] ). The 
practical applications of such expansions, however, have been most exploited by 
molecular physicists and chemists, where the multipole moments of molecules with 
various symmetries have been measured for decades (for example [2401 157]). and by 
geophysicists, with models of the magnetic fields of the planets within the Solar System 
readily found in the literature (for example |2bTj ). 

In this section we derive the magnetic field components of an arbitrary multipole 
of order I in spherical tensor form directly from the magnetostatic potential (where 
I — 1, 2, 3, 4, 5, . . . represent the dipole, the quadrupole, the octupole, the hexadecapole, 
the dotriacontapole, and so on). We show how the field components can be expressed 
in terms of either the polar, or the equatorial, field strength of the particular multipole 
being considered. We extend the analysis to demonstrate how the field components can 
be modified via the inclusion of an additional boundary condition designed to mimic the 
effects of plasma opening field lines to form a stellar wind, and show how the field can 
be written in co-ordinate free form. For stellar or planetary applications it is natural to 
employ spherical tensors, and later choose spherical polar components of the magnetic 
field vector B. However, in light of recently published work by the authors of [160J 
and [161] . who follow |148j by adopting a Cartesian tensor approach for the multipole 
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expansion, we conclude this section by demonstrating how an alternative approach to 
the multipole expansion can lead to differing expressions for the large scale magnetic 
field components. 

Throughout this work we assume a standard spherical polar coordinate system 
(r, 6, 0) with the coordinate origin (r = 0) taken to be the centre of the star (or 
equivalently a planet), < 9 < ix and < <fi < 2tt, with 8 = corresponding to the 
stellar rotation pole. When Cartesian coordinates are considered, the stellar rotation 
axis is assumed to be aligned with the z-axis and = the rc-axis. 



3.1. Magnetostatic expansion 

We are interested in deriving expressions for the multipolar field components that can be 
used to describe the large scale magnetospheres of stars (or equivalently planets). Their 
external magnetic fields are generated due to the dynamo action and the distribution 
of current sources internal to the star/planet. Expansions of the magnetostatic scalar 
potential, which we denote for a finite continuous current distribution can be carried 
out in several ways - by introducing Debye potentials [B2], using spinors |246| . or most 
elegantly, directly from Maxwell's equations for a magnetostatic field [83] (see also [2l]). 
In cgs units the magnetostatic Maxwell equations are, 

V-B =0 (3.1) 

47T 

V x B = —J, (3.2) 

c 

where J is the source current density, assumed localised near the coordinate origin 
in figure [I] By taking the curl of both sides of (3.2) and using a vector identity for 



V x (V x B) and (3.1) it is straightforward to show that 



Air 

V 2 B = -— VxJ. (3.3) 

c 

In a region which is source-free (for example in the region external to the star/planet) 
the complete form of the field B can be determined purely from its radial component 
r • B (see appendix B of [82], and [841 . for general proofs). By operating on both sides 



of (3.3) with r- and using another vector identity and (3.1) it can be shown that, 

47T 

V 2 (r-B) = r-VxJ. (3.4) 

c 

This is Poisson's equation, which has a solution in terms of the static Green's function 
for the Laplacian, |r — r'| _1 , 

r .B = l/r>; v ' x ,, J V 

c J |r — r | 

where r denotes a field point external to the star, r' denotes a source point internal to 
the star, while J' = J(r') and the source J is also assumed to be to internal to the star, 
as illustrated in figure [TJ 
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Figure 1. A current source J internal to a star which also contains the origin of the 
coordinate system 0. r denotes a distant field point external to the star, where the 
potential due to the current source at source point r' is to be calculated. 



By applying the cosine rule to the triangle in figure [I] the |r — r'| 1 term in (3.5) 
can be re-written as 

-1/2 



1 + 



2 ( - | cos6 
r 



(3.6) 



where 6 is the angle between r and r'. The term in the square brackets is the generating 



function for Legendre polynomials, which allows (3.6) to be be re-written as 

Ir-rr^E^costf). (3.7) 
l 

Using the addition theorem for spherical harmonics, which expresses the Legendre 
polynomials Pi(cos6) as the sum of the product of the spherical harmonics Yj 



and Y* m {9, 0) over the range m = —I, . . . , I, equation (3.5) can be rewritten as, 

r • B = I f drV • V x J'^E (^Tl) r^ Ylm(9 '> ^ (3 ' 8) 

The spherical harmonics are given by, 



Y, 



Im \ 



2/ + l\ 1/2 f(l-m)\\ 1/2 



An 



Pi m (cos8)e 



im<f> 



(3.9) 



(/ + m)\ 

for m > [while for m < 0, Y l{ _ m) {0 , (j>) = (-l) m Y* m {6,<f))]. The definition of the 
spherical harmonics differs between research areas via the inclusion or omission of the 
first two bracketed terms, or parts thereof. The definition we use here includes the 
Condon-Shortley phase [the (— l) m term]. Other definitions lead to differing expressions 



for the associated Legendre functions P/ m (cos#), defined below by (3.21). Defining the 
non-primitive spherical magnetostatic multipole moments asjjj 

j 1 • '~' Ur < al *y-v x j', (3.io) 



Mi 



i m 



c(l + 1) 



dr'r n Y, 



Im \ 



The difference between primitive and non-primitive multipole moments is discussed in Appendix B 
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allows (3.8) to be written as 
r B 



EE 

l m 



+ 



47T 



21 + 1 



M lm Y* m {9A)/r 



i+i 



(3.11) 



where the reason for the inclusion of the additional (I + 1) factor in (3.11) will become 



obvious during the integration of the separable differential equation (3.14) below. We 



note that he non-primitive multipole moments, defined by (3.10) can be re- written in 



several equivalent ways using various vector identities [82]. Alternatively, the |r — r 



'i-i 



term in (3.5) can be expanded in a Taylor series and then written in traceless form, 



in analogy with the electrostatic Cartesian tensor derivation in Appendix B Interested 
readers can find details of this in [85] and [86J. 

External to the star, in the source free region, (3.2) further reduces to V x B = 0. 



This condition can be satisfied by writing the field B in terms of a magnetostatic scalar 
potential 

B = -V*. (3.12) 

In order to determine the field components B r and Bg, required to describe the large 
scale structure of a stellar (or equivalently a planetary) magnetosphere, we first need to 



derive an expression for Operating on both sides of (3.12) with r- gives, 



r B 



Or 



(3.13) 



By equating (3.13) and (3.11) a separable partial differential equation is created that 



can be solved for 

dm 

dr 



EEc 



47T 



21 + 1 



M lm Y?JB,<j>)lr 



i+i 



(3.14) 



With the assumption that the potential ^ — > as r — > oo, (3.14) can be integrated to 



give, 



EE 



Air 



21 + 1 



M lm Y* m {9A)/r 



i+i 



(3.15) 



where the integration constant is zero (see [85J for further discussion on this subtle 
point). Equation (3.15) gives the general form of the multipole expansion of the 



magnetostatic potential in spherical tensor form. Note that the correct number of 
components (i.e., 21 + 1) of the non-primitive multipole moment of order I, M\ m 
with m = —I, . . . , +1, occurs automatically using spherical tensors (compare this with 



the Cartesian tensor method briefly discussed in £ Appendix B). The corresponding 



expression for the spherical tensor form of the non-primitive electric multipole moment 
Qi m is obtained from M im in (3.10) by replacing r' • V x J'/[c(l + 1)] with the charge 
density p(r'). 

In this paper we are interested in deriving expressions for the axial multipole field 
components, which correspond to the m = terms of (3.15) when we choose space-fixed 



axes with z along the multipole moment symmetry axis. The potential then cannot 
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depend on the azimuthal angle 0, so that only terms with m = can contribute to 



(3.15). By substituting for the spherical harmonics using (3.9) the potential becomes 



4tt \ 1/2 

^( r ) = ^i^rn) M lQ P l {co & 9)/r l+l . (3.16) 



The Miq terms are determined directly from (3.10) 

1/2 

c{l 



m 10 = -y^—j J dv'r' 1 {^pj Pticosey • V x J' 



where the quantities Mi are defined as the magnetic multipole moments, with Mi = 
the dipole moment, Mi = Q the quadrupole moment, M3 = f2 the octupole moment, and 
so on. The Ith component of the scalar potential of the large scale stellar magnetosphere 
is therefore given by, 

*, = ^*Kcos0). (3.18) 

We note that in general M/ m is a complicated function of the orientation of the 
current source distribution, but for axial distributions it is a simple function of the 
orientation (8, <p) of the symmetry axis [H], i.e., 

M lm = M l Y lm (6,<P). (3.19) 

A corresponding expression can be derived for the Cartesian form of the non-primitive 
multipole moment [871189] . and analogous spherical and Cartesian tensor formulae exist 
for the non-primitive electric multipole moments [87, 88J. 

3.2. The magnetic field components 

As the stellar magnetic field is described by the gradient of the magnetostatic scalar 
potential, B = — V\P, the field components in spherical coordinates of a multipole of 
order / are obtained via, 

Br = ~^ * = (3.20) 

while for the axial multipoles that we consider in this paper = 0. Figure [2] illustrates 
a field vector B at a point along a field line, decomposed into the B r and Bg components. 
Noting that the associated Legendre functions can be written as 

P lm (^) = (l-x 2 r /2 ^Pi(x) (3.21) 
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Figure 2. The left hand panel shows a field vector B decomposed into the radial B r 
and polar Bg components at a point along a field line a distance r from the centre 
of the star, at a co-latitude of 6. The field components are used to illustrate their 
definitions and are not to scale. The right hand panel shows the unit vectors in the 
radial and polar directions and M; is a unit vector along the symmetry axis of axial 
multipole 



where x = cos# and Pi(x) are the Legendre polynomial^] 

1 d' 



P,(x) 



2 l l\ dx l 



(3.22) 



the magnetic field components of an axial multipole of order I can be obtained from 

Mi 



(|3~20} using (13181) and flOIL 



B r 



+ 



,1+2 



MtPA cos 9) B e 



„l+2 



-Pa (cos 



(3.23) 

Pi(cos8) and Mi is the multipole 



where we have made use of the fact that Pi (cos9) 
moment. 

The field components for the lower order multipoles can then be obtained. For the 
dipole magnetic field, 



B r 



2/i 



cos 6 Ba 



sin 9. 



The axial quadrupole magnetic field components are given by, 
3Q 



B r 



2r 4 



(3 cos 2 6-1) B e 



— cos sin 0, 



and the axial octupole field components are, 
2VL 



B r 



(5 cos 2 9 - 3) cos 6 B e 



3Q 
2r5 



(5cos 2 #- l)sin0. 



(3.24) 



(3.25) 



(3.26) 



Expressions for higher order multipoles can be derived from (3.23). 

For some applications, including lunar magnetism, the multipolar B fields within 
the source region are required, which can be conveniently represented in terms of contact 

% We remind readers that the Condon-Shortley phase, a term of the form (—1)™, has already been 
included in our definition of the spherical harmonics, and is therefore not included as a pre-factor in the 
associated Legendre functions. Users of IDL should note that the legendre package of IDL does include 
the (— l) m term. It is also worth noting that the geophysics community uses Schmidt quasi-normalised 
associated Legendre functions, by convention, which leads to slightly different expressions for B r and 
Bg, see [262] . No such convention exists in the astrophysics community. 
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fields involving the Dirac delta function 8(r) [89]. Because V x B / in the source 
region we cannot use the usual magnetic scalar potential \I> in this region. Instead B can 
be represented by the two scalar Debye potentials ip and \ [BS], B = L^ + Vx L%, where 
L = — ir x V is the angular momentum operator. The two terms are the toroidal and 
poloidal components, respectively, of B. The coefficients of contact field terms for ip, 
X and B involve the primitive magnetic multipole moments. An application of contact 
multipolar fields to lunar magnetism is described in reference [HE]- The corresponding 
electrostatic multipolar contact fields |88j , which can be derived using the electric scalar 
potential since V x E = in the source regions, are mentioned in §Appendix B| 



3.3. General expressions for magneto spheres 

3.3.1. Polar field strength Rather than specifying the strength of the various multipole 
moments, it is more convenient to discuss the polar strength of each component, B l f ole , 
i.e. the strength of the particular multipole component at the stellar rotation pole. At 
the rotation pole, r = i?* and 8 = 0, P/(cos#) = 1 for all values of /, and the field of all 
the axial multipoles is purely radial. From (3.23) any multipole moment can be written 
as 

Mi = R l + 2 B l f ole /(l + 1) (3.27) 

where B l f ole is the polar field strength of the /th order multipole. The dipole, quadrupole 
and octupole moments can be written as \i = Bl> pole R^/2, Q = B 2 ' pole R^/3 and 
Q = B^ pole R^/4: respectively, which allows the field components to be re-expressed 
in a more convenient form, 

B r , m =Bl' pole (^j cosO, (3.28) 

B e , dip =±Bl^(^j sine (3.29) 

B rmad = l -Bl*° le (3cos 2 fl-l), (3.30) 

B e , quad = B 2 ^ (^j cos e sin e, (3.31) 

B r , oct = l -Bl' pole (^j (5 cos 2 £-3) cos (3.32) 

B e ,oct =\Bl' pole (^j (5cos 2 #-l)sin#, (3.33) 

For an /th order multipole the general field components can be derived from (3.23) using 
( p7) , 

(D \ l + 2 T)l,pole / D \ l + 2 

^1 Pi(cose) B e = j^ T (^\ P n (cose). (3.34) 
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These simple expressions for the magnetic field components of axial stellar (or planetary) 
magnetospheres are straightforward to adapt as inputs to numerical simulations. They 
can be used to derive the components of an individual Ith order multipole, while 
linear combinations of the various multipoles may be used to construct expressions 



for more complex multipolar fields. The use of (3.34) requires knowledge of Pi(cos6) 
and Pn(cos8) as well as the polar field strength of each multipole component being 
considered. The Legendre polynomials and associated Legendre functions can be looked 



up in tables, or derived through a combination of (3.21) and (3.22) and use of Bonnet's 



recursion formula. For models of stellar magnetospheres the polar field strength of 
each multipole component is determined observationally by decomposing the field into 
poloidal and toroidal components, each of which is then expressed as a spherical 
harmonic expansion. The coefficients of such a fit to the data contain information 
on the strength of the individual field components [4"9] . 



3.3.2. Equatorial field strength It is also possible to derive expressions for B T and 
Bq which include the stellar equatorial, rather than the polar, field strength. Various 
authors define the strength of a low order multipole as being the field strength at the 
stellar equator (for example [TUl |3l] for a dipole and |173] for a quadrupole) while others 
follow the convention used in the previous subsection, and define the strength as the 
polar value (for example |175[ 194"]). With 9 = n/2 and r = it can be seen from (3.28) 



and (3.29) that the equatorial field strength of a dipole is 1/2 of the polar value. For a 
quadrupole, it can be shown from ( |3.30[ ) and fl3.3l| ) that B%' equ = B} pole /2^\ However, 
this is not a general result for higher order multipoles. For example, using the expressions 



for the octupole field components, (3.32) and (3.33), we find that B^ equ = 3B^ pole /8. In 



Appendix A we derive a general relationship between Bl ,equ and B l f ole for a multipole 



of arbitrary /. 

At the stellar rotation pole the field is purely radial, and we therefore end up with 



single expressions for B r and Bq and the particular multipole moment Mi [see (3.27) and 



(3.34)] that are valid for all multipoles when expressed in terms of B l f ole irrespective of 
the I value. This is not the case if the expressions for B r , Bq and Mi are re-expressed 
in terms of the equatorial field strength B[ ,equ . For odd I number axial multipoles the 
field in the stellar equatorial plane has only a polar (6) component (B r = 0,Bq ^ 0), 
while for even I number axial multipoles the field only has a radial component in the 
equatorial plane (B r ^ 0, Bq = 0). Therefore the expressions for B r , Bq and Mi, when 
written in terms of the equatorial field strength, are different for odd and even / number 
multipoles. For odd I number multipoles, such as the dipole, the octupole and the 



dotriacontapole, our expressions (3.34) can be re-written in terms of the equatorial field 



strength B\ equ rather than the polar field strength B l f ole using the results developed in 

+ Note that for the quadrupole, in the equatorial plane, Bg — and B r — — B 2,pole /2. This then gives 
an equatorial held strength of B^ equ = {B 2 + B 2 ) 1 ' 2 = B^ pole /2. 
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Appendix A, 



Br = MiiOT^ (B.\ fl(cM „ (3 . 35) 



l\ * \r , 

B, = 2-[('-l)/2]![( i+ l)/2] !i? ^ CM' +2 p , l(cos() ), (3.36) 



(Z + l)! * W 

while for even / number multipoles, such as the quadrupole and the hexadecapole, the 
corresponding expressions are, 

*>l\(1 /OM12 / d \ ; + 2 



Aiim 2 B ie qu ( R*\ Pl{cos6) (3.37) 
l\ \ r J 



mm 2 1+2 



Be = -f^h B * eQU {v ) P ^ cos ^ ( 3 - 38 ) 

The constant terms can be modified by introducing double factorials, however, this 
does not lead to a significant simplification. Clearly the expressions for B r and Bg when 
written in terms of the equatorial field strength, lack the simplicity and elegance of the 



corresponding expressions based on the polar field strength (3.34). 



3. 3. 3. Incorporating open field by using a source surface boundary condition References 
[2] and [219J introduced the source surface boundary condition in order to produce global 
potential (i.e. current free) field extrapolations of the Sun's coronal magnetic field 
from maps of the photospheric field. This outer boundary condition of the potential 
field source surface (PFSS) model mimics the effect of the solar wind dragging and 
distorting the field lines of the solar corona, and gives a simple way to incorporate open 
field into global magnetospheric models. In reality the distortion of the field by the 
coronal plasma will induce a current system, and therefore a proper solution for the 
field structure requires a solution to the equations of MHD. None-the-less PFSS models 
have been used extensively in the study of solar magnetism. At some height above the 
solar surface, the source surface Rs, the plasma pressure in the corona pulls open the 
field lines forming a wind. Above the source surface the field is purely radial and is 
often described by a Parker spiral |189j . The source surface represents the radius at 
which all of the field becomes radial, but there is no reason why it may not do so closer 
to the surface. For the Sun Rs is typically taken to be 2.5 R , a value consistent with 
satellite observations of the interplanetary magnetic field (see the discussion in |212j ). 
This same boundary condition (with differing values of Rs) has since been applied to 
stellar field extrapolations for young rapid rotators [110] and pre-main sequence stars 
|91j . providing a simple method of incorporating open field regions along which a wind 
could be launched. In this section we show how our general expressions for B r and 



Bg, (3.34), are modified for a stellar magnetosphere with a source surface (note that 
PFSS models have also been applied to planetary magnetospheres, for example, |223j ). 
In §[4] we compare the results of PFSS extrapolation models with more complex MHD 
simulations. 
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Figure 3. Field lines for the lowest order multipoles (a dipole, I = 1, a quadrupole, 
I = 2, and an octupole, I — 3) with a source surface set to Rs ~ 3.4i?* plotted as 
the dashed line. The field line shape is calculated using the components B r and Be 
obtained from (3.48) and (3.49). For higher order multipoles there are regions of open 
field at lower latitudes along which a stellar wind could be launched. The magnetic 
fields are symmetric about the x and z axes, while a multipole of order I has 21 shells 
of closed loops around the entire star (for example 



The large scale magnetospheric field must satisfy Maxwell's equation that the field 
be divergence free, (3.1), and this, combined with (3.12) (the PFSS model assumes 
that B is source free), means that the magnetic scalar potential must satisfy Laplace's 
equation, 

V 2 * = 0. (3.39) 
This has a separable solution in spherical coordinates of the form 

* = E + bi m r~ {l+l) ] P lm (oos 6y m * (3.40) 



where here we change from the normalised spherical harmonics in (3.9) to conform 
with our previously published models of stellar magnetosphere |110[ I9~T| I9~4] , and where 
the coefficients a\ m and h\ m are determined from the boundary conditions. The first 
boundary condition is to specify the radial component at the stellar surface. For 
field extrapolation models, this is determined directly from the observationally derived 
magnetic surface maps. In our case, as we are considering axial multipoles, the radial 
field at the stellar surface (r = R*) is given by (3.34), 

B r {R«) = B l f ole P l {cosd), (3.41) 

for a multipole of order /. The second boundary condition is that at the source surface 
Rs the field becomes purely radial, 

Be(R s ) = B^Rs) = 0. (3.42) 



The magnetic field components themselves can be derived from (3.12) using (3.40) 
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From (3.44) and (3.45) it is clear that boundary condition (3.42) is satisfied if 
bim 
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(3.46) 



while for the axial multipoles 



m 



0) it can be seen from (3.45) that B 



■ 0. From 

(3.40) it can be seen that (3.46) is equivalent to the assumption that \l/(r = R$) is 
an equipotential surface. Substituting ( 3.46[ ) into (3.43), with m = 0, and applying 



boundary condition (3.41) and noting that Pi (cos9) = Pi( c os9) } an expression for c% 



in terms of R* and Rs can be derived for a particular multipole of order /, 

TDl,pole 

am = * 7tt^- (3-47) 

Substituting this result and ( 3.46[ ) into (3.43) and ( |3.44[ ) with m = gives, for a 
particular / value, general expressions for B r and Be (valid in the region R* < r < Rs) 
for the large scale magnetosphere with a source surface, 
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(3.49) 

where in deriving (3.49) we have used the fact that Pu(cos9) = —dPi(cos9)/d9 [see 
1]. The field lines of the lowest order multipoles, with a source surface, 

rd9/B e 



(3.21 ) with m 



are sketched in figure pj These are found by numerical integration of dr/B r 



with B r and Bg given by (3.48) and (3.49), although an analytic solution can also be 



found, but we do not discuss this here. Note that the magnetic field components with 



the imposed source surface boundary condition are the same as (3.34) multiplied by 



correction terms. Changing the radius of the source surface Rs modifies the structure 
of the entire magnetosphere, with more open field, along which a wind could be launched, 
available for smaller values of Rs- Equations (3.34) are recovered in the limit of 
Rs — > oo. 



3.3.4- Coordinate free field components and tilted magneto spheres The initial ZDI 
results on V2129 Oph and BP Tau have suggested that the octupole field component 
of the magnetospheres of accreting T Tauri stars contains a significant fraction of the 
magnetic energy [501 El]- The dipole components of their magnetospheres, however, 
remain the most dominant at typical disc truncation radii, despite the large scale dipole- 
like field being distorted close to the surface of the star [91]. Composite magnetic fields 
consisting dipole plus octupole field components have been used for many years by the 
solar physics community in the study of coronal mass ejections (for example (6j,EQj). For 
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stellar magnetism, [162] and |211j have recently presented MHD simulations of accretion 
to stars with composite dipole-octupole fields, the latter comparing their model directly 
with observations of V2129 Oph |50j. In their 3D models the octupole and dipole 
moment symmetry axes are tilted relative to each other, and to the stellar rotation axis, 
and the three axes do not lie in one plane. Their prescription for the total field B' of 
axial multipole I is presented in coordinate free form. The total field B' can be written 
using our equations ( 3.23[ ) as 



B' 



J+2 



(/ + l)P,(cos0)f + P n (cos6)d 



(3.50) 



Let M; be a unit vector along the symmetry axis of axial multipole I. It can be seen from 
figure [2] (right panel) that M; • f = cos 9 and Mj ■ = — sin 6*. To make the derivation 
easier to follow we have assumed in figure [2] that is aligned with the stellar rotation 
axis, and that both lie in the same stellar meridional plane (for example, the a^-plane). 
However, the results derived in this section apply generally to tilted multipole symmetry 
axes, and with appropriate coordinate and vector frame transformations are equally 
applicable to the case of two (or more) axial moments with arbitrary tilts with respect 
to the stellar rotation axis. As M; can be written generally as = cos^f — sin 9Q 
then (Mj • 9)9 = — (Mj ■ ?)r. Thus a simple expression for 9 can be derived, 
9 = — cosec^M; + cot#f. Using this result to eliminate in (3.50) then the total field 
may be written as 
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(3.51) 



In Appendix A the associated Legendre functions (with m = 1) and the Legendre 
polynomials are written as series expansions. Using (A.l) and (A.2) to replace Pi(cos6) 
and Pn(cos6) in (3.51), and using (3.27) to replace the Z-th order multipole moment 
with the polar field strength, then after some manipulation, a general result for the 
total field of an arbitrary titled axial multipole of order I in coordinate free form can be 
obtained 
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where N = 1/2 or N = (I — l)/2, whichever is an integer. This general result can be 
used to construct models of composite magnetic fields consisting of the fields of two or 
more axial multipoles, arbitrarily tilted with respect to the stellar rotation axis, such 
as presented by |162j . Analogous to (3.52) a similar expression can be derived for a 
magnetosphere with regions of open field introduced by applying the source surface 
boundary condition. Starting from (3.48) and (3.49) and following the same argument 
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used in deriving (3.52) we obtain, 
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3.4- Difference between a spherical and Cartesian tensor approach 



The authors of [1601 1161] have recently constructed models of composite dipole- 
quadrupole magnetic fields, however, the expressions that they derive for the quadrupole 
field components are a factor of 1/2 smaller than those derived in this paper (see 



£3.2). In their work they follow the approach of |148j who develop an expression for 
the magnetostatic potential exterior to the star due to a pseudo magnetic "charge" 
distribution interior to the star. In their approach the magnetostatic potential is 
derived by analogy with the electrostatic potential expanded using Cartesian tensors. 



In Appendix B we consider the electrostatic case by expanding the potential for a 
finite static charge distribution in Cartesian coordinates. As part of that derivation 



the non-primitive quadrupole moment must be defined (see equation (B.9)). There 
are three different definitions of the traceless quadrupole moment tensor used in the 
literature (for example |139l 11481 1129] ). which ultimately leads to different expressions 
for B r and Bq, and explains why the expressions used in [160] and [161j are a factor 
of 1/2 smaller. However, as demonstrated in [87], the factor of 1/2 arises naturally 
when the potential is expanded using spherical harmonics. Given that stars and their 
circumstellar environments, and planets, are most straightforwardly described using a 



spherical coordinate system, (B.9) is the most convenient definition of the non-primitive 



quadrupole moment. Furthermore, as we now show, our equations (3.34) represent the 



simplest way of expressing high order field components, and do not suffer from any 
ambiguity that can arise due to the differing definitions of multipole moments in use. 

As pointed out in Appendix B three different definitions for the non-primitive 
(traceless) electrostatic quadrupole moment tensor Q are usec^*" 

Qi = l^Qii^iTi 



rfl) 



Q2 = ^( 3r i r i - rf 1 ) 

% 

i 

where is the tensor product of the vectors r-j, I is the second rank identity tensor, 
and Q : T^ 2 )(r) used below is the double dot product representing the full contraction 

Further definitions of Qi are possible and are discussed in [82 for both the electric and magnetic 
multipolcs in terms of the equivalent spherical tensor Qi m . In |86| equivalent Cartesian forms of the 
traceless magnetic quadrupole moment corresponding to Qi arc discussed. 
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of the tensors Q and the gradient tensor T^(r) = VV(l/r) (see Appendix B for 
full details). Following [148] . the authors of |160j 1161] adopt the Q 2 definition of the 
non-primitive quadrupole moment. As a result of this, it is trivial to show that the 
quadrupole potential, and therefore the field components B rma d and Bg tqua d derived in 
[1601 1161] , are a factor of 1/2 smaller than we have derived in this paper, based on the Qx 
definition. Adoption of the Q3 definition results in a third different form for the potential 
and field components for the quadrupole. However, it is easy to demonstrate that our 
general expressions (3.34), based on the polar field strength, produce quadrupole field 
components that are independent of how Q is defined. 

We define the non-primitive quadrupole moment for the electrostatic Cartesian 
tensor case between (B.7) and (B.8). The (electrostatic) quadrupole potential can 
clearly be written as 5>2( r ) = : T^ 2 ^(r), where the constant 77 depends on whether 
the Qi, Q2 or Q 3 definition is used (with Qi,r) = 1/3, with Q2,?7 = 1/6 and with 
Q 3 ,?7 = 1/2). Following the argument in Appendix B.2 and carrying out the double 
dot product to determine the electrostatic potential $2 (which has an identical form to 
the magnetostatic potential derived from equation (3.18)), and then using (3.20) to 
determine B r and Ba, we find that, 



B r = ^t1](3 cos 2 9 — 1) Bq = —^?7 cos 9 sin 9. 



(3.54) 



2 r 4 • r 
At the stellar rotation pole, where the field only has a B r component, it is clear that 
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As both B r and Bq are directly proportional to both 77 and Q, and as Q itself is inversely 
proportional to 77, then the field components are independent of 77 and consequently 
independent of the chosen definition of the non-primitive quadrupole moment tensor. 
Thus, our equations (3.34) represent the simplest way of expressing the field components, 
and do not suffer from the ambiguities encountered in the literature due to differing 
definitions of the various multipole moments. 



4. Magnetospheric accretion models with multipolar magnetic fields 

Models of the magnetospheres of accreting T Tauri stars have traditionally assumed 
dipolar magnetic fields. In the past few years however, motivated in part by the 
availability of stellar magnetic surface maps, new models have been developed which 
consider fields with an observed degree of complexity. Two types of models which 
incorporate multipolar stellar magnetic fields have been developed. Firstly, those which 
assume that the large scale magnetosphere can be modelled as a potential field, and 
then include effects of the coronal plasma, and secondly, 3D magnetohydrodynamic 
simulations. 
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4-1. Development of PFSS models and comparison with MHD field extrapolations 

Over the past 40 years the PFSS model has undergone several modifications. Reference 
[197] considered how the coronal gas pressure distorted the large scale field. In order to 
account for the effects of the coronal plasma, the authors of [222] introduced an analytic 
field extrapolation model with a non-spherical source surface, but only considered the 
dipole component of the solar field. The resulting source surface was a prolate spheroid 
with major axis aligned with the solar rotation axis. Reference |153j subsequently 
extended the ideas of Schulz to incorporate more generalised fields (see also |224j ). It 
is interesting to note, however, that the prolate spheroid source surface also arises in 
more complex MHD field extrapolation models, at least at solar minimum |212j . At 
solar maximum, when the photospheric field is more complex, the source surface is, on 
average, more spherical. 

Various other studies have examined the effects of including current sheets and 
other volume currents into solar magnetic field extrapolation models. Such models 
better reproduce the topologies of coronal streamers detected during solar eclipse 
observations, as well as satellite observations of the large scale heliospheric field (for 
example \220\ I263[ El 12701 1271[ ITT]). As such detailed observations are not available for 
forming solar-like stars we do not discuss the various models here. Interested readers can 
find up-to-date discussions comparing PFSS models with the results of more complex 
current sheet source surface (CSSS) models in [225] and [116j . 

MHD field extrapolation models are more desirable as they not only include 
information on the magnetic field structure, but also about how fields interact with, and 
are influenced by, outflows and the coronal plasma. For the solar corona, the first global 
MHD simulations to incorporate magnetic maps of the photospheric field were developed 
in the late 1990s (for example [250] and [157] ). A thorough and balanced overview of the 
advantages and disadvantages of the use of the PFSS model compared to MHD models 
is presented in [212] (see figure [4]). The PFSS allows 3D global field extrapolations to be 
produced quickly with moderate computing resources. In contrast, MHD models require 
powerful computing resources. The main advantage of MHD models is their ability to 
include the non-potentiality of the field and time dependent surface transport effects 
(meridional circulation, differential rotation, supergranular flows etc). PFSS models 
produce static field configurations and cannot incorporate time dependent effects, such 
as magnetic reconnection events. It remains an open question, however, how important 
surface effects are on T Tauri stars, although the authors of [57] have recently reported 
the detection of surface differential rotation on the low mass CTTS V2247 Oph. 

MHD models require the inclusion of a number of free parameters, particularly 
with regards the thermodynamic properties of the coronal plasma. In the MHD model 
of [212] a simple polytropic equation of state for the coronal plasma is considered, 
with the temperature and density at the base of the solar surface treated as free 
parameters. Other authors have considered more complex implementations of the energy 
equation; reference [58] for example, considers heating and cooling terms in the MHD 




Figure 4. A comparison of the field structure of the solar corona obtained using the 
potential field source surface method (PFSS; right hand panel) with a MHD solution 
(left hand panel) at solar minimum (top panel) and close to solar maximum (bottom 
panel). For both the PFSS and MHD models the field is extrapolated using the same 
starting positions at the photosphere with field lines assigned the same (arbitrary) 
colour to aid comparison. Reproduced by permission of the AAS, and Dr. P. Riley 
from EH]. 



energy equation. Unfortunately, such models also introduce more unknown physical 
parameters, the choice of which directly influences the resulting field topology. In 
contrast to MHD models, the PFSS model produces a unique solution for the coronal 
field |3j. Although reference |212] concludes that the PFSS model often produces results 
close to the more physically realistic MHD models, when considering the large scale 
coronal structure, the authors anticipated that discrepancies between the models would 
increase once new solar vector magnetograms became available (readers are referred to 
|213] and |244] for some of the latest work on solar field extrapolations). Extrapolations 
of the local magnetic field in active regions have recently been compared to the observed 
3D field structure derived from the STEREO satellite data. Reference |218j concludes 
that the potential field approximation is poor at reproducing the observed fields, while 
reference |158j concludes that provided the active region under consideration is relatively 
stable then the potential field approximation is adequate. 

It is important to remember however that small scale magnetic features, such as 
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bipolar groups, that are easily detected on the Sun, remain well below the achievable 
resolution of stellar magnetic maps. A star in the Taurus star forming region, at a 
distance of ~ 140 pc [128J, is some 29 million times further from the Earth than the 
Sun. Stellar spectropolarimetric measurements can only probe the large scale properties 
of the magnetosphere, and not the small scale field regions where the PFSS model most 
frequently breaks down; thus the PFSS model has since been extensively applied to 
produce coronal field extrapolations from stellar magnetograms derived through ZDI 
(for CTTS, however, there is additional non-potentiality induced due to the interaction 
of the large scale magnetosphere with the disc, as we discuss in £4.3). 

In addition to the finite achievable resolution Zeeman-Doppler maps of stellar 
surface fields are also subject to missing information due to inclination effects, whereby 
for a given stellar inclination, much of one hemisphere is hidden to an observer. All of 
these limitations, and their effects on the derived coronal structure and X-ray emission 
properties, have been thoroughly examined in [125J. Of the three limitations (finite 
resolution, the suppression of the Zeeman signal in dark spots, and the unobservable 
hemisphere) they find that the largest effect on the global coronal structure is caused by 
stellar inclination effects. What theoretical models assume for the flux distribution in 
the hidden surface area changes the way the larger scale field lines connect to opposite 
polarity regions between the hemispheres. However, the authors of |125] conclude 
that the field in the visible hemisphere is reliably reproduced by field extrapolation 
models, and in particular, the quantities calculated from PFSS models of T Tauri 
magnetospheres (the location of hot spots, accretion filling factors, disc truncation radii 
etc) are unaffected by missing information. For CTTS the detection of high latitude 
accretion spots strongly suggests field configurations that are antisymmetric with respect 
to the centre of the star. If it were the case that the field in the hidden hemisphere 
was such that stellar magnetosphere was symmetric with respect to the centre of the 
star, high latitude accretion spots would be difficult to explain. The authors of |125] 
also conclude, however, that missing flux in magnetic surface maps causes PFSS models 
to overestimate the amount of open flux. This may have important implications for 
T Tauri angular momentum loss (see ^ in the following section we discuss the PFSS 
model applied specifically to T Tauri coronae). 

A problem for PFSS models of stellar coronae is where to locate the source surface. 
Unlike solar observations, in-situ satellite measurements of stellar heliospheric fields 
are not available. Indirect indicators of coronal structure have therefore been used to 
estimate the maximum extent of the closed field regions of stellar coronae. For AB Dor, 
a star for which several magnetic maps have been derived [45], 147], slingshot prominences 
have been detected extending to ~ 4.5R*, in excess of the equatorial corotation radius 
of ~ 2.7R* [36, 47J. Thus, in the first PFSS field extrapolation models to be applied to 
a stellar magnetic map (of AB Dor), |110j made a conservative estimate of Rs = 3.4R*. 
However, changing the location of the source surface influences the structure of the entire 
corona, in particular the extent of open field relative to closed field regions. Using the 
AB Dor magnetic maps [lllj extended the PFSS model to include the effects of coronal 
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plasma. Assuming that stellar coronae are isothermal and that the plasma trapped 
along the closed loops is in hydrostatic equilibrium, enables the gas density at each 
point within the corona to be determined. 

In the Jardine coronal model if the gas pressure exceeds the magnetic pressure at any 
point along a field line loop, that field line is assumed to be torn open, and the coronal 
plasma lost in a wind. Such loops would therefore be dark in X-rays. In such a way 
the model allows the stellar X-ray emission properties such as the global X-ray emission 
measure, or the amount of rotational modulation of X-ray emission, to be estimated 
using the coronal fields extrapolated from magnetic surface maps. With a large source 
surface, however, often the simulated values of the coronal density are too low, suggesting 
that more compact coronae, with small values of Rs are required. At first, this appears 
to contradict the requirement of large extended magnetic structures for prominence 
support. However, even if the closed field regions are confined close to the surface of the 
star, extended stable magnetic loops may form due to the reconnection of wind-bearing 
open field lines |112j . This model provides a mechanism for supporting prominences 
beyond corotation, and beyond any reasonable estimate of the source surface radius. 
The location of the source surface in stellar models remains a somewhat free parameter, 
although changing its value affects the amount of modulation of X-ray emission and the 
X-ray luminosity predicted. The coronal model does, however, include an additional free 
parameter. It is assumed that the gas pressure at the field line foot points (equivalently 
the density) scales with the magnetic pressure (p oc Bq). Adjusting the free parameter, 
the constant of proportionality between these two pressures, affects the total amount of 
X-ray emission predicted. For smaller values the gas pressure is reduced at each point 
along the coronal loops and therefore the magnetic field is able to contain more of the 
coronal plasma. For larger values, the opposite is the case, and more loops are unable to 
contain the plasma and are assumed to be blown open. For T Tauri stars specifically the 
authors of |113j developed a method of constraining the free parameter by comparing 
the model predicted X-ray emission measures with those derived from a large sample of 
stars in the Orion Nebula Cluster, for each assumed surface magnetic field distribution. 
We note that such free parameters are not unique to stellar coronal models, and are 
commonly employed for model of solar loops, for example [260J. 

Independently obtained X-ray observations of stellar coronae can provide a test of 
the coronal fields derived via field extrapolation. In order to test the model properly, 
however, contemporaneous spectropolarimetric (from which the magnetic surface maps 
are derived) and X-ray satellite observations are required, see [106J. The simultaneity 
of the multiwavelength observations is crucial as if there is a large separation between 
the ground based spectropolarimetry and the space based X-ray observations there is 
a danger that the magnetic field of the star will have evolved significantly during the 
delay. Thus, the X-ray properties predicted from the derived magnetic map may be 
significantly different from those observed. This of course would not provide a true 
test of the theoretical models, as it would be impossible to ascertain if the difference 
was due to magnetic evolution of the stellar corona, or inadequacies of the model, or 
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both. However, the lack of significant change in the large scale field topology of the 
accreting T Tauri stars BP Tau and AA Tau apparent from spectropolarimetric data 
taken years apart, suggests that a strictly simultaneous observing strategy, although 
highly desirable, may not be required [5U [52] . A contemporaneous observing strategy 
was successfully employed in |1U6] on AB Dor and is currently being used in order to test 
the ability of PFSS models to capture the true magnetospheric geometry of accreting T 
Tauri stars [95J, which we now discuss. 

4-2. Potential field models of T Tauri magneto spheres with complex fields 

The PFSS model was extended through a series of papers in order to construct models 
of the magnetospheres of accreting T Tauri stars [901 EH EH EU 1113] . The initial 
simulations used magnetic fields extrapolated from surface magnetograms of young 
rapidly rotating zero-age main sequence stars, AB Dor and LQ Hya, as pre-main 
sequence maps were not available at the time [901 I9~T] . By adjusting the stellar 
parameters to typical T Tauri values the extrapolated fields were divided into three 
distinct regions. Close to the star the complex and loopy field lines contained the X-ray 
emitting corona, while regions of open field (typically at high latitude) were available 
to carry a stellar wind. The larger scale field, that which interacts with the disc, is 
simpler in structure and more "dipole-like" . This simple larger scale field is, however, 
distorted close to the star by the complex surface field regions [91]. This model for 
T Tauri magnetospheres assumed that any of the large scale field lines that passed 
through the stellar equatorial plane would support gas accretion from the disc on to 
the star, provided they passed through the disc interior to the corotation radius down 
to some inner disc truncation radius. Interior to corotation the effective gravity points 
inwards towards the star and accretion would naturally take place. This simple model, 
the first to incorporate multipolar magnetic fields into models of the accretion process 
on to T Tauri stars, successfully reproduced the observed correlation between X-ray 
luminosity and stellar mass [198, 1113] . the observed rotational modulation of X-ray 
emission [691 192] , the observed reduction in X-ray luminosity of accreting stars relative 
to non-accretors [H3], and the correlation between disc mass accretion rate and stellar 
mass [91]. However, it is worth noting that despite several publications describing the 
possible origin of the accretion rate - stellar mass correlation (originally discovered by, 
although thus-far not credited to, reference [102] ). it may be nothing more than an 
artefact of observational selection and detection limitations [35J. 

Figure [5j shows potential field extrapolations from photospheric maps of the 
magnetic fields of the accreting T Tauri stars V2129 Oph and BP Tau. There is a clear 
distinction between the simple large scale field, and the more complex surface field. 
From the rotational modulation of unpolarised and circularly polarised profiles of the 
Call infrared triplet emission (and in particular the contribution from accretion spots 
at the foot points of accretion flows), |50j and [51] infer that the majority of accreting 
gas flows into spots at high latitude (see [51] for details of the construction process for 
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Figure 5. Numerical field extrapolations of the magnctospheres of the accreting T 
Tauri stars V2129 Oph and BP Tau (left and right panels respectively). Open field 
lines are coloured blue, with closed fields line in white. 

excess Call emission maps, which has recently been improved and further developed 
by the authors of |52j). In order to ensure that the large scale field is able to reach 
high latitudes, the source surface in such models must be set to at least the equatorial 
corotation radius R co . This does not mean, however, that the disc need be truncated 
at corotation. By assuming that circumstellar discs are truncated where the differential 
magnetic torque due the stellar magnetosphere is balanced by the differential viscous 
torque in the disc, reference [HI] calculated disc truncation radii using the extrapolated 
fields of BP Tau and V2129 Oph. For BP Tau, a completely convective star with strong 
dipole and octupole components, the disc was found to be truncated at ~ 0.7R co . 
The disc of V2129 Oph, a star which despite its young age has already developed a 
small radiative core and which has a dipole component four times weaker than that 
of BP Tau, the truncation radius was ~ 0.5-R co . The process of accretion of gas from 
well within corotation should exert a spin-up torque upon the star in the absence of an 
efficient angular momentum removal mechanism. We note, however, that based on a re- 
evaluation of the accretion related emission lines in the optical spectrum of V2129 Oph 
its mass accretion rate may be an order of magnitude lower than previously considered 
by both [5TI] and [HI]. With a lower accretion rate the torque balance calculation would 
result in a larger disc truncation radius. A stronger dipole component and/or a weaker 
octupole component would have the same effect. 

The incorporation of complex multipolar magnetic fields into magnetospheric 
accretion models naturally provides an explanation for small accretion filling factors 
(the fractional surface area of the star covered in accretion hot spots). Magnetospheric 
accretion with an aligned dipole field tends to produce filling factors that are typically 
an order of magnitude larger than observationally inferred (for example [25 2] ) . For 
magnetospheric accretion with non-dipolar magnetic fields accreting gas is funnelled on 
to discrete regions of the stellar surface [HDl HO]- The prediction of small accretion 
hot spots spanning a range of latitudes is not unique to PFSS models of T Tauri 
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magnetospheres. Hot spots at low latitudes are a natural consequence of considering 
non-dipolar, as well as tilted dipolar, magnetospheres. They also arise in MHD models 
of magnetospheric accretion (as we discuss in the following section) and in a recent 
generalisation of the Shu X-wind model |179j . 

The multipolar X-wind model of reference [179] is constructed from 2D 
axisymmetric potential stellar magnetic fields. Only the odd Z-number multipoles, which 
have B r = in the equatorial plane, are considered^] As with the dipolar X-wind model, 
all field lines initially originate from the starjff] During the initial phases of evolution 
it is assumed that matter somehow opens the large scale closed field threading the disc. 
This produces a disconnected region of open field at the corotation radius from where an 
outflow is launched. Although a non-dipolar magnetosphere results in funnel flows from 
the disc arriving at the star at different latitudes, the basic properties of the X-wind 
model remain unchanged, namely the assumption of disc-locking and that of trapped 
flux within a small region (the X-region) at the corotation radius. The field lines both 
interior and exterior to corotation are pinched together, while angular momentum from 
the funnel flow is transferred backwards through the X-region and removed by the X- 
wind. Thus, the multipolar X-wind model, as well as the dipole version |2291I187] , allows 
accretion to take place without a net transfer of angular momentum to the star. The 
process of accretion proceeds without a spin-up torque being exerted on the star, which 
would have slowly increased the stellar rotation rate towards break-up speed. This is 
desirable as most accreting T Tauri stars have rotation rates about an order of magnitude 
below the break-up rate (for example p2]). The model is admirable in that it provides a 
simultaneous description of both accretion and outflows, processes which are intimately 
linked in accreting T Tauri stars (large collimated jets, for example, are not observed in 
non-accreting systems). However, the multipolar X-wind model suffers from the same 
limitations as its dipolar counterpart (see Furthermore, the multipolar X-wind 

model, nor the PFSS field extrapolation models, incorporate time dependent effects, 
such as the field evolution due to the star-disc interaction, or magnetic reconnection 
events. The modelling of such complex processes requires a full MHD solution. 

4-3. 3D MHD models of T Tauri magnetospheres with non- dipolar fields 

Reference [21 lj presents 3D MHD simulations of V2129 Oph. This represents the 
culmination (thus far) of a series of papers by the same group of authors examining the 
star-disc interaction with ever increasing physical reality - from 2D MHD simulations 
[206J, to 3D models with tilted dipole magnetospheres |207[ 1208] (see also |136j ). 
simulations with composite dipole-quadrupole magnetic fields |16U] , then with composite 
magnetic fields where the dipole and quadrupole moments are tilted by different amounts 
relative to the stellar rotation axis, and in different planes |161j . and finally similar 3D 

jj Reference |179) erroneously and repeatedly refers to the quadrupole as being an odd multipole and 
the I = 3 and / = 5 multipoles as the quadrupole and the octupole. In fact, such Z-numbers represent 
the octupole (I = 3) and the dotriacontapole (1 = 5). 

ff The model has recently been extended via the inclusion of a separate disc magnetic field |141j . 
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Figure 6. A MHD simulation of the magnetic field of the CTTS V2129 Oph showing 
the distortion of an initially potential field (t — 0, left hand panel) due to the star-disc 
interaction. The middle and right panels show two different orientations of the field 
structure after 14.4 days (t = 8 in simulation units), a little more than two stellar 
rotation periods. The upper panel is the field closest to the star, which shows little 
departure from the initial potential configuration. The lower panel shows the formation 
of a magnetic tower, which arises due to the differential rotation of the field line foot 
points anchored in the disc and on the stellar surface, and being wrapped around the 
rotation axis. Figure reproduced, with permission, from |211j . 



models but with composite fields consisting of a dipole plus an octupole field component 

The most recent simulations have shown that the potential field approximation 
remains valid in regions where the magnetic stresses dominate over the material stresses 
of the accreting gas |162tl211j . as was argued previously |lf 5j . The 3D MHD models can, 
however, include the time dependent effects of the star-disc interaction. The larger scale 
field is found to quickly depart from a potential field configuration due to the shearing 
of the field lines caused by the differential rotation between where they are anchored on 
the stellar surface and where they thread the disc, see figure [6j which leads to inflation 
of the field. The field lines wrap around the rotation axis forming a magnetic tower 
( |211j . and references therein). 

The shape and distribution of accretion hot spots at the base of flows of accreting 
gas are particularly sensitive to the field topology. The inclusion of multipole moments 
of different strength, tilt relative to the rotation axis and the other moments, as well 
as the plane of the tilt, are all important in determining the location of hot spots. For 
example, in the case of accretion to dipole-octupole composite fields, a strong dipole 
component dominates the kinematics of the in-falling gas, leading to the formation of two 
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discrete funnel flows for a slightly tilted dipole moment [162] . In other simulations with 
a stronger octupole component, the more complex field regions dominate leading to the 
formation of lower latitude spots. Intriguingly, the position and shapes of photospheric 
hot spots determined from the simulations appear to remain fairly constant in time 
[161] . This immediately suggests that the commonly observed complex variability of 
accreting T Tauri stars, evident from their photometric light curves, is due to unsteady 
(clumpy) accretion flows, with a variable mass accretion rate. Alternatively this complex 
variability may be due to the restructuring of the stellar magnetic field due to dynamo 
processes interior to the star (see the discussion in [162j )|f] 

The authors of [16U] and [161] have also calculated the angular momentum flux 
on to the star due to accretion and the outward angular momentum flux along (polar) 
stellar open field lines. They find that the main spin-up torque on the star is due to 
the magnetic connection with regions of the disc interior to corotation, with a smaller 
spin-down torque arising from the outward angular momentum flow along the open field. 
The spin-up torque from the accreting matter itself is a factor of 20 — 100 times smaller 
than the spin-up torque from the magnetic connection to the inner disc for a purely 
dipolar stellar magnetic field, but only 10 — 20 times smaller for a purely quadrupolar 
field [160] . In all cases the MHD models presented in references [160] and [161] feature 
stars that are experiencing a net spin-up torque. For stars with magnetic fields that 
are more complex than a dipole, the spin-up torque is less than for the dipole models 
despite the disc being truncated closer to the star. This may initially seem contradictory, 
however the authors of [161] argue that it arises due to the smaller connectivity with 
the inner disc than they find for dipolar fields. Thus the authors of [161J conclude that 
the spin-up torque experienced by a star with a non-dipole magnetic field may not be 
as severe as for stars with dipole magnetic fields. Further investigations incorporating 
generalised multipolar stellar fields, simulated for many rotations of the star would be 
welcomed to explore these suggestions further. 

The MHD models of accretion to non-dipolar magnetic fields do not yet incorporate 
restructuring of the stellar magnetosphere due to surface transport effects, or time 
variable accretion flows (but see [210] and [146] ). However, it is not yet clear 
observationally whether or not the surface fields of accreting T Tauri stars vary 
significantly with time. Reference [HI] derived two magnetic maps of the accreting 
T Tauri star BP Tau from circularly polarised spectra taken approximately 10 months 
apart. This corresponds to almost 39 stellar rotations, and while the maps appeared to 
show an apparent quarter phase shift in the field (most likely related to a small error in 
the assumed rotation period [21]), the large scale field structure is remarkably similar at 
both epochs. The same result was found for AA Tau where magnetic maps derived from 
spectropolarimetric data taken about one year apart also show little difference in the 
large scale field topology [52] . This would then suggest that time variable mass accretion 
has a larger role to play in T Tauri photometric variability; however, the authors of 

f It is also worth noting that the authors of [161] conclude that stars which show simple sinusoidal 
light curves can also host complex non-dipolar magnetic fields. 
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[186], who search for variations in accretion rates over short (~ hours - days) and long 
(~ months) timescales conclude that accretion hot spots rotating across the star is the 
dominant cause of variability. This result also appears to be consistent with earlier 
line profile variability studies of individual stars that suggested that accretion flows are 
modulated with the stellar rotation (for example, [77 | 11171 |119[ [T6"]). However, it should 
also be noted that the common detection of large X-ray flares believed to extend to 
several stellar radii may indicate large scale field evolution [721 E3]- This, coupled with 
the recent detection of strong differential rotation on the completely convective T Tauri 
star V2247 Oph [57] , suggests that more spectropolarimetric data on a greater number 
of stars is required to thoroughly investigate the issue of evolving magnetic fields. 

An important finding from the MHD simulations is that accretion hot spots are 
inhomogeneous, being hotter and denser in their cores compared to their peripheries, as 
most of the kinetic energy of the accreting gas is carried in the centre of the funnel flow (a 
finding first discussed for accretion along tilted dipole magnetospheres - |208] ). Thus, hot 
spots appear smaller when viewed at higher temperature and densities (i.e. at shorter 
wavelengths). This has important implications as it suggests that UV observations used 
to estimate accretion filling factors under estimate the true size of the spots (being 
sensitive only to the hotter regions). This also appears to explain why accretion filling 
factors and mass accretion rates derived from density sensitive line triplets detected in 
high resolution X-ray spectra are found to be an order of magnitude smaller than those 
derived from optical/UV spectra, with the X-ray emission arising from only the hottest 
regions of the spot [981 E] • However, it must also be remembered that unsteady accretion 
flows may alter the accretion filling factor, [7], for example, find a factor of more than 40 
difference in the accretion filling factor for BP Tau. This is of course a variability study 
of a single star, and long term observations targeting the same stars repeatedly would 
be welcome to search for variations in accretion spots sizes and distributions. The non- 
uniform variability of accreting T Tauri stars clearly presents a formidable challenge 
for theoretical models (see the discussion in |122j ). Future models that incorporate 
magnetic surface maps as a boundary condition, and which consider stellar surface 
transports effects, will drive further progress in this field. 

5. Summary and applications to outstanding problems 

Over the past few years instrumentation has advanced to the stage where the magnetic 
fields of stars across the Hertzsprung-Russell diagram can be probed in unprecedented 
detail. Of particular interest are the magnetic fields of forming solar-like stars, as they 
allow us to study the history of the Sun at the epoch of the formation of the Solar 
System. The ability to produce surface maps of their magnetic fields is a remarkable 
achievement, given that stars (with the exception of the Sun) are mere points of light 
in the night sky. However, those points of light are rotating, and by monitoring the 
distortion of Zeeman signatures as stars rotate, ZDI studies have revealed the complex 
nature of T Tauri magnetic fields. In this work, we have demonstrated how results 
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developed in the classical electromagnetism and molecular physics literature can be used 
to derive analytic expressions for the axial multipole (m = 0) magnetic field components 
of a stellar, or equivalently a planetary, magnetosphere. The resulting expressions (3.34) 
depend only on the polar field strength of the particular multipole component and on the 
associated Legendre function P fl (cos0) and Legendre polynomial p(cos#). Our general 
expressions for B r and Be in terms of the polar field strength are valid regardless of the 
definition of the non-primitive (traceless) multipole moments. 

Complex magnetic fields can be created by considering linear combinations of the 
multipole field components, for example, by following a similar approach to that used 
by the authors of [149J I131[ HblH 11611 1138] . who have all considered a dipole plus a 
quadrupole composite field, or the approach of the authors of |179] and |162j . who have 
considered higher order multipoles. Our expressions for B T and Bg provide the most 
straightforward method of constructing complex magnetospheric geometries consisting 
of axial multipoles. 

The analytic descriptions of axial multipoles derived in this paper, when used as 
inputs to new models, will increase our understanding of many important astrophysical 
processes. Applications include models of the magnetic star-disc interaction, wind 
launching on both the pre-main sequence, and main sequence, and studies of the 
magnetic interaction between stars and orbiting close-in giant planets. 

The interaction between T Tauri magnetospheres and the disc is believed to control 
the rotational evolution of the star [20] and may also be responsible for the collimation 
and launching of outflowing winds and jets [66]. This process may determine the mass 
accretion rate, setting the lifetime of the disc, and consequently the brief window of time 
in which planets may form. The star-disc interaction depends on the stellar magnetic 
field topology, with the disc truncation radius being particularly sensitive to the relative 
polar strength of each of the multipole field components [94J. Recent observations 
of the classical T Tauri star V2129 Oph revealed a complex magnetic topology with 
field modes up to I = 15 required to fit the data [50]. The magnetic energy was 
concentrated dominantly in the octupole field component, with a weak dipole component 
dominating at larger radii [94J. As an alternative to numerical field extrapolations (see 
[115] and [51]) the large scale magnetosphere of V2129 Oph may be (rather crudely) 
approximated analytically by a combination of the / = 1 and I = 3 multipoles, giving 
B r = B rdip + B TOCt and B 9 = B e dip + B e octl where "dip" and "oct" refer to the dipole 
and octupole components respectively. Each of the terms in such linear combinations 
are easily derived from (3.34). Such composite expressions for the radial and polar field 
components do not violate Maxwell's equation that the field be solenoidal (V ■ B = 0). 

MHD simulations of the star-disc interaction which incorporate observationally 
derived magnetic surface maps will represent the next major advancement in the field. 
However, this will require significant computational resources. The magnetic surface 
maps contain information about many high order field components, thus any simulations 
will require small grid resolutions at the stellar surface to handle the steep gradients in 
field strength, and this must be coupled with a large enough computational domain to 
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capture the behaviour at the disc interaction region. It is questionable whether or not 
the resulting simulations will yield new insights not already gained from simulations with 
tilted dipole plus octupole composite fields |162[ 1211] . Such composite fields broadly 
match the magnetic topologies obtained to date through ZDI of accreting T Tauri stars 
[501 [M]. A more fruitful line of research may be the inclusion of surface transport effects, 
differential rotation, supergranular diffusion and meridional circulation, into MHD 
simulations, as well as a more complete consideration of accretion shocks, and surface 
waves, generated by the high velocity impact of the dense accreting gas |137] l38| 139] . 

Another anticipated application for our general expressions for B r and Bg, is 
in models of stellar winds. Understanding the rotational evolution of stars at all 
evolutionary phases requires knowledge of how the stellar wind torque varies with the 
stellar parameters (for example |lU9j ). Recently, the authors of [174] have considered 
how the stellar wind torque for main sequence solar-like stars depends on quantities such 
as the stellar radius, mass, and rotation rate, the mass outflow rate, and the equatorial 
field strength of the stellar magnetosphere. By considering both dipolar and quadrupolar 
magnetospheres, they find that the stellar wind torque is particularly sensitive to the 
assumed field geometry. However, as pointed out [174] , further simulations are required 
to fully quantify the effects of varying the field topology. Our analytic expressions for 
the magnetic field components provide an easy way of incorporating more complex field 
geometries in to stellar wind models. 

Exoplanet migration may be influenced by the stellar wind plasma, as well as by 
how the stellar magnetosphere interacts with the disc. Simulations [209] . and analytic 
work |156] l70j. suggest that the inner disc hole, cleared by the star-disc interaction, may 
provide a natural barrier that decreases the rate of inward migration of forming planets. 
However, if a multipolar magnetosphere were to be considered, the structure of columns 
of accreting gas within the magnetospheric gap may alter the migration rate of planets 
[209J. By considering the azimuthal ram pressure of a stellar wind, the authors of [164J 
find that a planet will migrate inwards or outwards depending on whether the stellar 
rotation period is greater or less than the planets orbital period. The current picture is 
confused, however, as the authors of [254] have presented 3D MHD simulations of stellar 
winds and conclude that they have little influence on the migration of exoplanets. But, 
as pointed out [164] . multipolar stellar magnetic fields are expected to have a non- 
uniform distribution of open field, along which a stellar wind could be launched. A 
migrating planet, in the stellar equatorial plane, may therefore be immersed directly in 
the stellar outflow. This potentially will have a larger influence upon planet migration 
than would be expected from dipole stellar field models where the wind is launched from 
the star in a direction away from the equatorial plane. Such arguments, however, are 
highly speculative and more quantitative work in this area is required. 

Other anticipated applications for our analytic expressions for multipolar 
magnetic fields include models of magnetic star-planet interaction, and of planetary 
magnetospheres. The magnetic fields of Jupiter and Saturn, for example, are well 
known to be multipolar with estimates of their multipole moments readily found in the 
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literature (for example [261]). Outside the Solar System, it has been argued that orbiting 
close-in giant exoplanets can lead to stellar activity enhancements. This is an emerging 
research area, which may provide a method for characterising exoplanetary magnetic 
fields, and consequently their internal structure [228J . Variations in the activity of host 
stars, synchronised with the planetary orbital period, have been observed, although 
sometimes with a short phase lag, perhaps due to the planet perturbing the stellar field 
lines (for example [258]). A possible explanation for this is that a planet ploughing 
through a stellar magnetosphere triggers a release of energy stored-up in the coronal 
field by decreasing its relative helicity (for example |150[ I151j ). However, such models 
have yet to consider multipolar stellar magnetic fields. Realistic multipolar fields have 
significantly less magnetic flux at the position of orbiting close-in planets, even if the 
large scale topology is "dipole-like" in structure [H] . More observations and theoretical 
modelling are required to confirm if magnetic star-planet interaction can indeed lead to 
stellar activity enhancements. A more promising method of investigating exoplanetary 
magnetic fields is through the detection of planetary radio emission, arising due to the 
electron-cyclotron maser instability, where electrons are accelerated along the field lines 
of the planetary magnetosphere due to the electric field generated from the reconnection 
of the stellar and planetary magnetic fields |115j . 

The current generation of spectropolarimeters are providing the community with 
unrivalled new datasets with which to probe stellar magnetism as a function of stellar 
age and spectral type. Such datasets can provide crucial guidance for new theoretical 
models designed to deepen and broaden our understanding of all of the above physical 
processes. In turn, such new models may be used to make predictions that can be tested 
with future instrumentation. The HARPS spectropolarimeter will provide better access 
to astrophysical objects in the southern hemisphere sky [237] . and, as with the PEPSI 
(Potsdam Echelle Polarimetric and Spectroscopic Instrument) at the Large Binocular 
Telescope, will allow Stokes Q and U linear polarisation studies, important for the 
investigation of magnetic fields at the inner edge of circumstellar discs [2411 1108] . In the 
longer term, SPIRou, a nIR spectropolarimeter for the Canada-France-Hawai'i telescope 
to be implemented in 2014, will provide simultaneous Zeeman broadening measurements 
and magnetic topology information. 
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Appendix A. Relations between the equatorial and polar field strength for 
a multipole of arbitrary order / 



The Legendre polynomials are defined by (3.22), however, an alternative definition is 
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where A" = 1/2 ot N = (I — l)/2, whichever is an integer [142] . Using this alternative 
definition expressions for B r and Bg based on the stellar equatorial, rather than the 
polar, field strength, can be derived. 

Appendix A.l. Odd- order multipoles 

For multipoles of odd-order /, for example, a dipole (I = 1) or an octupole (I = 3), 
AT = (Z — l)/2. In the stellar equatorial plane B r = 0, and the field only has a Bg 



component. From (3.23) the Bg component depends on the associated Legendre function 



Pn(x). Calculating Pn(x) requires the Legendre polynomials to be differentiated [see 



(3.21)]. From (A.l) and (3.21) we obtain 
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Therefore the field strength at the stellar equator for an odd I number multipole is, 
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By noting that the only possible negative term in (A.4) is the (—1)^ term, and 
that |(— l)( z-1 )/ 2 | = 1 regardless of the odd I value, (A.4) can be re-arranged for the 



multipole moment Mi. Comparing the result to (3.27), a relationship between the polar 



and equatorial field strengths can be determined for an odd /-number multipole, 
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These result is not valid for even /-number multipoles. 
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Appendix A. 2. Even-order multipoles 

For multipoles of even-order I, for example, a quadrupole (Z 
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1/2 and (A.l) can be re- written as 
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where we have substituted for N 
equator of the star, where 6 = tt/2, the only non-zero term of (A.6) is the final term. 
For even / number multipoles the field is purely radial in the equatorial plane (Bg = 0). 
Therefore, from (3.23), 
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The equatorial field strength for an even I number multipole is then, 
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Noting that the only possible negative term in (A.8) is the (— l)'/ 2 term, and that 
|(— l) z / 2 | = 1 regardless of the even I value, (A.8) can be re-arranged for the multipole 
moment Mi. By then comparing the result to (3.27), a relationship between the polar 
and equatorial field strengths for an even /-number multipole can be determined, 
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This expression is not valid for odd /-number multipoles. 



Appendix B. Electrostatic expansion using Cartesian tensors 

The authors of |160] have recently presented a derivation of the quadrupole component 
of the potential expansion in Cartesian coordinates, based on the electrostatic approach 
presented in |148] , but using pseudo magnetic "charges" . In this Appendix we derive an 
expression for the electrostatic potential $ external to a volume (for example, the star) 
containing the "charges" and the coordinate system origin. As part of the derivation, for 
terms in the potential expansion corresponding to the quadrupole and the higher order 
multipoles, a choice must be made regarding the definition of the traceless multipole 
moment. Various different definitions are used in the literature, which ultimately leads 



to different expressions for the magnetic field components (see the discussion in £3.4). 

We consider a finite discrete charge distribution with N charges g« contained within 
a volume which also contains the coordinate origin. The electrostatic potential $ at a 
distant field point r due to a charges at source points r 1; r 2 , . . . is given by 

*W = £iFrFr (R1) 

i 1 1 

This, and the other results below, are straightforward to adapt for a continuous charge 
distribution by replacing the sum over the individual charges with an integral over 
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volume of the charge density (i.e. g» — > J p(r)dr). The multipole expansion in terms 
of Cartesian tensors can then be obtained by expanding |r — r^j" 1 as a Taylor series in 
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where = 1/r, T^(r) = V(l/r), T^(r) = VV(l/r), ... are gradient tensors. The 
primitive multipole moments relative to the origin are given by q = £\ q%, M- = Q'i 1 "?, 
= ^i<7i r i r i> O = Sj% r i r « r i an d so on 5 where q is the total charge, and \i, 8 and 
O are the dipole, quadrupole and octupole primitive moments respectively. The term 
primitive moment [199] is used here to distinguish these multipole moments from the 
traceless multipole moments (which are referred to in the literature as the multipole 
moments) discussed below in ^Appendix B.2 The Cartesian tensors which define the 
primitive moments are symmetric (e.g. Q a p = 9p a ) but not traceless (e.g. E a QQ ^ 0). 
In this notation quantities such as represents the tensor product of the vectors 
(in this case specifically the dyadic product). Note that a tensor rr has components 
(rr) a/ 3 = r a Tp and similarly (rrr) Qj g 7 = r a r^r 7 etc. The double and triple dot products 

are represented by : and : respectively and represent full contractions of the appropriate 
rank tensorsjf For example, : T^ 2 ^(r) = ^2 a ^2pd a pT^ represents a full contraction, 
yielding of the second-rank tensors and T^ 2 ). (Since the tensors used in 

this paper are all symmetric we can also write : T^(r) = d a pT^).) Equation 



(B.3) is valid for an electrostatic potential generated by true electric charges [57], or for 
a magnetostatic potential generated by pseudo magnetic charges. 



The multipole expansion (B.3) is valid for field points r outside the source region. 

^88] the multipole contributions to the electric field inside the 



For some applications 
source region can be expressed in terms of so-called contact terms, involving the delta 
function S(r). It is interesting to note that the primitive multipole moments are 
necessary to fully describe these contact field terms; the traceless moments are sufficient 
for the external (long-range) multipole fields (as we discuss below for the electric case 



and in {3.1 for the magnetic case) but not for the internal (contact) multipole fields. 
Magnetostatic multipolar contact fields [S9] are mentioned in 33.21 



1 r I 



while a full contraction of two 



| The tensor product is also denoted in the literature by 

rank-n tensors A and B can also be written as A-n-B or A[n]B. For example, O : T/ 3 ' = O'3-T™ = 
0[3]T( 3 ). Partial contractions of tensor products can be defined similarly [57] but will not be needed 
here. 
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Appendix B.l. The dipole term 



The dipole term of the multipole expansion, the second term in (B.3), is straight-forward 
to derive. The first gradient tensor can be written as T^^(r) = V(l/r) = — fr~ 2 where 
f is a unit vector along the direction of r. The dipole moment vector (x has three 
components, however, by transforming to the principal axes frame such that the z-axis 
lies along the direction of yi, gives a dipole moment of \i = fi z (principal axes frames are 
discussed in greater detail below for the quadrupole moment). By performing the dot 



product contained in the second term of (B.3) we obtain, 



$1 



cos 9. 



(B.4) 



where we have used the fact that z/r = cos 9, and where <3>i is the dipole component of 
the potential expansion. 



Appendix B.2. The quadrupole term 

As a symmetric second rank tensor, the primitive Cartesian quadrupole moment 9 has 
six independent components. However, the independent components can be reduced to 
five as terms of the form AI can be added to 9, where A is an arbitrary scalar and I the 
second rank identity tensor, since 



AI:T (2) (r) = AV 2 ^ =0 (r>0), 



(B.5) 



i.e. AI contributes nothing to the field external to the distribution. Note that in general, 
a primitive multipole moment of order I, which is symmetric, has (/ + 1)(Z + 2)/2 
independent components, and 1(1 — l)/2 traces |176j . Generally the non-primitive 
moment has (2/ + 1) independent components, which emerge automatically in the 



spherical tensor method (see the discussion of equation (3.10) in £3.1). The form of 
A is chosen to be (see for example |152j ). 



A 



(B.6) 



from which it is straight forward to show that 9 + AI is traceless, thus reducing the 



number of independent components to five. The quadrupole term in (B.3) is then 



* 2 (r) = \e : T«(r) = \ J> (r«r, - : T< 3 >(r) 



iQ:T( 2 )( 



where the traceless (or non-primitive) quadrupole moment is defined as, 
Q = l^Qii^iTi -rfl). 



(B.7) 
(B.8) 

(B.9) 



The derivation of the corresponding octupolar term $3(1") in terms of the non-primitive 
or traceless Cartesian octupole moment is discussed in reference [87] and in the arXiv 
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version of this paper. At this stage it is worth pointing out that there are three 
different definitions of the non-primitive quadrupole moment used in the literature 
(for example [1391 I148[ 1129] ). which ultimately leads to different expressions for B r 
and Bg (see £3.4). However, the definition we chose here ensures that the Cartesian 
tensor approach produces expressions for the electrostatic case, and equivalently the 
magnetostatic potential case (see below) expanded using a spherical tensor approach, 
that are entirely equivalent. The components of the traceless quadrupole moment tensor 
are written as, 

Q a P = o ^ft( 3r «* r i/9 ~ r » 2( W)> (B.10) 

i 

where 8 a p the Kronecker delta (6 a p = 1 if a = (3 or =0 if a ^ (3). 

Multipole moments are mathematical constructions, the values of which, in general, 
depend upon the choice of origin |176] 1199] . The first non- vanishing multipole moment 
is independent of the origin [87]. For a body of general shape three of the quadrupole 
moment components represent the orientation of the body-fixed axes (of the star or 
planet) relative to the space-fixed axes of the coordinate system being used to define 
the components [87]. As the multipole moment tensors are real and symmetric it is 
always possible to transform to the so-called principal axes frame, which has the effect 
of reducing all off-diagonal elements of ( |B.9 ) to zero. This is equivalent to calculating the 
principal moments of inertia (moments of a mass distribution) in mechanics problems 
(see [79] for a thorough discussion of the principal axis transformation). Choosing the 
body-fixed axes to coincide with the principal axes reduces the number of independent 
components of Q to two (recall that as Q is traceless, Q xx + Q yy + Q zz = 0). In the 
axisymmetric case the number of independent principal axes components reduces to 
one. To see this choose the body-fixed z-axis of the star (the rotation axis) as the 
symmetry axis so that Q xx = Q yy by symmetry and both are equal to —Q zz /2 since Q 
is traceless. For axial distributions the body-fixed component Q zz = Q is referred to as 
the quadrupole moment. 

In order to derive the final expression for the quadrupole term of the potential 
expansion, we also require the components of the second rank gradient tensor. Noting 
that T( 2 )(r) = VV(l/r) = (3ff — I)r~ 3 (r ^ 0) then the components of T^ 2 ^(r) can be 
written as, 

T$ = (3r a r p -r 2 8 a p)r- 5 . (B.ll) 



Calculating the double dot product (where Q : = ^2 a YlaQafiTg2) in (B.8) gives 
the quadrupole component of the multipole expansion for the axial case, 

$2 = |^(3 cos 2 #-1). (B.12) 

The expressions above for the dipole and quadrupole terms of the electrostatic 
potential expansion derived from a Cartesian tensor approach, $i and $2 as given by 
equations (B.4) and (B.12), are entirely analogous to the equivalent terms obtained 
form the spherical tensor approach to the magnetostatic expansion, and ^2 as 
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derived from equation (3.18). Derivation of the traceless octupole and higher order 
traceless multipole terms of the potential expansion using the Cartesian tensor approach 
is more cumbersome and time consuming; however, the components of the non-primitive 
octupole and hexadecapole moments (1 = 3 and I = 4 respectively) can be found in the 
literature (for example [240] ) , and reference [S7| also provides a general expression for the 
l-th traceless moment of such an expansion. Electrostatic multipole expansions can also 
be carried out using spherical tensors by expanding the |r — r'| _1 term in equation (B.l ) 
and using the addition theorem for spherical harmonics following a similar argument to 
that discussed in §3.1| The resulting expression for <E> Z , the general term of electrostatic 
expansion, is entirely analogous to the general term of the magnetostatic expansion 
given by equation (3.18). A detailed first principles derivation can be found in [87J. 



Appendix C. Octupole term of the scalar potential 



(This appendix is included in the astro-ph version only and will not appear 
in the published version of the article). 

The octupole moment tensor O, as a third rank tensor, has three possible traces 
O a 0p). In general, a primitive multipole moment of order I, which is symmetric, has 
(I + 1)(Z + 2)/2 independent components, and 1(1 — l)/2 traces |176j . The non-primitive 
moment has (21 + 1) independent components in general [as seen automatically using 
spherical tensors - see (D.4)], or seven for the octupole moment (I = 3) as we show 



below. Following a similar argument to that of the quadrupole case, see Appendix 



B.2, a term of the form (al)^" 1 , where a is an arbitrary vector and (al) sym denotes a 
symmetrised tensor with components (al)^™ = (a a 8^ + + a y 5 a p), can be added 
to the primitive tensor since 

(aI)«w»: T (3)( r ) =Qj ( r >0). (C.l) 
Because a is arbitrary we can add a sum of such terms of the form 

i 

such that all three traces of the new octupole moment tensor vanish. The octupole term 
in (B.3 ) is then, 

$ 3 (r) 



-ioiT^r): 
o 

15 V ; 



r^r, - l -r} (r,I)^ m 



: T (3) ( 



where f2, the traceless octupole moment tensor, is defined as, 

i 

This has components, 



n 



a/37 



p(r)dv(5r a r f3 r^ - r 2 [r a S^ + r p 5 ia + r 7 <5 Q/3 ]) 



(C.3) 
(C.4) 

(C.5) 
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where we have switched to considering a continuous charge distribution rather than a 
summation over individual charges (Y^ p(r)dr where p(r) is the charge density) 



for ease of notation. We have chosen the numerical prefactors in (C.3) and (C.4) so that 



(C.4) agrees with the definition using spherical coordinates. For an arbitrary charge 



distribution (C.5 ) has seven independent components in a general reference frame. In the 



axisymmetric case, where there is one independent principal axes component, called the 
octupole moment Q = Vt zzz , there are seven non- vanishing principal axes components, 



i.e. 



n 



yij-- 



zyy 



with Q 2 



The components of the third rank gradient tensor, T^ 3 )(r) 
I)r _3 l, can be written as, 



VT< 2 )(r) 



V[(3ff 



T 



(3) 



a/37 



(-3r 7 )(5r Q r /3 r 7 - r 2 (r a 5^ + rp8 ia + r 7 <5 a/3 )). 



(C.6) 



Thus carrying out the full contraction with the traceless octupole moment tensor 
(fi : T^ 3 ) = J2 a S/3 ^atpiT^a) gives an expression for the octupole term of the 



multipole expansion in the axial case 

n 



2r 4 



(5 cos 2 6-3) cos 9. 



a 



(C.7) 

is equal to Q3 defined by 



The octupole moment introduced in this section Q 

J rfrp(r)r 3 P 3 (cos 9). The axial magnetic octupolar field has the same 
M 3 where M 3 is obtained from the definition of Q3 by replacing 



(D.6), i.e. Q 3 



form as (C.7), with f2 



p(r) by r • V x J/ (4c). 



Appendix D. Electrostatic expansion using spherical tensors 

(This appendix is included in the astro-ph version only and will not appear 
in the published version of the article). 

Derivation of the components of the non-primitive octupole and hexadecapole moments 
(1 = 3 and I = 4 respectively) can be found in the literature (for example |24U] ) . 
however, for the dotriacontapole (I = 5) moment and beyond, the Cartesian tensor 
approach quickly becomes cumbersome and time-consuming. Multipole expansions can 
also be carried out using spherical tensors, which is more straightforward. A detailed 
first principles derivation can be found in |87j . 



By applying the cosine rule to the triangle in figure Dl the potential $(r) can be 
re-written as 

r 9 1 -V 2 

q . 



$ r 



E 



E 



1 + 



cos 9i 



(D.l) 



where 9 t is the angle between r and r,. The second term in the square brackets is the 



generating function for Legendre polynomials, which allows (3.6) to be be re-written as 

J 



i I 



(D.2) 
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Figure Dl. Coordinates and notation for a finite charge distribution inside a volume 
which contains the origin of the coordinate system 0. r is a distant field point where 
the potential due to discrete charges qi at source points is to be calculated. Only 
the first charge q\ is shown. 



Using the addition theorem for spherical harmonics, which expresses the Legendre 
polynomials Pi(cos9) as the sum of the product of the spherical harmonics Yi m (9i,(f>i) 



and Y, 



over the range m = —I, we obtain, 



$(V) 



EE 

l m 



An 



21 + 1 



Q lm Y* rn {e,4>)/r 



i+i 



(D.3) 



with Q\ m being the mth component of the spherical multipole moment tensor of order 
I 



Qim = y^g»r-y /m ( 



(D.4) 



where the spherical harmonics Yi 



lm\ 



are given by (3.9). The rotationally invariant 



form (scalar product) ^2 m Ai m Bf m of two spherical tensors occurs in (D.3), as it must 



and rotational invariance arguments can be used to give an elegant alternative derivation 
of (D.3) [87]. Note also that the spherical multipole moments correspond to the non- 
primitive (traceless) Cartesian moments for every I. Thus for / = 2, for example, there 
are five quadrupole components Qim, with m = —2, —1, 0, 1, 2. These can be written as 
linear combinations of the five independent Cartesian components Q a p 



For axial multipoles, as considered in this paper, we choose the space-fixed axes with 
the z-axis along the symmetry axis. The potential then cannot depend on the azimuthal 
angle 0, so that only terms with m = can contribute to (D.3), which therefore reduces 
to, 



where Qio is given by 



Q 



10 



21 



4:71 



1/2 



21 



47T 



1/2 



Qi- 



(D.5) 



(D.6) 



The quantities Qi are referred to as the multipole moments, with Q = q, Qi = /i, 
Q 2 = Q and Q 3 = Q with q the charge, and /i, Q and fl respectively the dipole, 
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quadrupole and octupole moments. The quantities q, \i and Q are defined in Appendix 



[B] and fl is defined in Appendix C In general Qi m is a complicated function of the 



orientation of the charge distribution, but for axial distributions it is a simple function 
of the orientation (6, 0) of the symmetry axis [87], i.e., 

Qi m = QiY lm (9,<f>). (D.7) 

An analogous expression for the traceless Cartesian multipole moments is given in 



The expressions (B.4), (B.12), and (C.7) derived using the Cartesian tensor method 



are obtained by using (D.5) and (D.6). The general term of the electrostatic potential 



expansion for the axial case is then, 



$i = ^Pi(cos6). (D.8) 
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